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CHAPTER 1 


INTRODUCTION 

The use of reduced-order models [1] and reduced-order observers [2] 
in the design of feedback controllers has been studied by several re- 
searchers. In addition, the development of eigenvalue/eigenvector as- 
signment techniques has received much attention in recent years. In 
this work, a reduced-order model is used with eigenvalue/eigenvector 
assignment techniques to design a constant state feedback controller for 
the original full-order system. The eigenvalues and eigenvectors con- 
tained in the reduced-order model are reassigned in the full-order sys- 
tem while those eigenvalues and associated eigenvectors not included in 
the reduced-order model remain unchanged in the full-order system. The 
constant state feedback matrix is implemented using output feedback with 
a reduced-order observer. It is shown that the eigenvalues and eigen- 
vectors of the closed loop full-order system remain unchanged when the 
reduced-order observer is implemented. 

1.1 Motivation 

During the past fifteen years significant advances have been made 
toward developing viable synthesis techniques for multivariable feedback 
control systems. Notable among these techniques is the eigenvalue/ 
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eigenvector assignment procedure. Early studies in this area focused on 
an algorithmic formulation of the spectral assignment by Srinathkumar 
[3], while later studies included a geometric formulation of the same 
problem by Moore [4], Kimura [5], and Davison and Wang [6]. Based on 
these theories, design procedures have been developed for approximating 
desired mode mixing [7], reducing eigensystem sensitivity to variations 
in plant parameters [8], reducing the effects of actuator noise on sys- 
tem performance [9] and modifying the resultant feedback gain matrix to 
specified gain constraints [10]. Recently, these procedures have been 
incorporated in a spectral assignment computer aided design package 
[11], A deficiency in all work concerning eigenvalue/eigenvector as- 
signment procedures is an absence of application of these techniques to 
real world design problems. A primary factor contributing to this prob- 
lem is the lack of understanding of how to use reduced-order models and 
reduced-order observers with spectral assignment procedures. 

Models representing the behavior of physical systems often consist 
of a very large nimber of coupled, linear differential equations. Such 
models are difficult to use when designing control systems due to exces- 
sive requirements for computer time and memory, and to the nimerical 
analysis problems inherently present when dealing with large systems of 
equations. It is, therefore, desirable to develop a design procedure 
which utilizes reduced-order system models. Simplification of large 
order dynamic systems has received the attention of many researchers in 
recent years. The major difficulty with this work is that only open- 
loop system behavior is approximated. Of concern when using reduced- 
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order models with eigenvalue/eigenvector assignment procedures is the 
fact that while the reduced-order model may approximate open-loop system 
behavior, the modeling error may be so great or of such a nature that 
actual closed-loop system performance is not acceptable. Also of con- 
cern is the closed-loop behavior of those modes of the original system 
which are not included in the reduced-order model. 

Full state feedback is implemented by the use of a dynamic observer 
system when there are fewer outputs than states. Since some states are 
usually available for measurement at the output, a reduced-order observ- 
er is desirable in order to minimize the complexity of the control sys- 
tem. Of concern is the effect of a reduced-order observer on the system 
eigenvalues and eigenvectors. 

1.2 Overview 

In this section an overview of the thesis is / given. A background 
of spectral assignment theory is discussed in Chapter 2. A subsequent 
design procedure implemented by Marefat in a computer aided design pack- 
age is presented next. This information provides a necessary foundation 
to support the material in the remaining chapters. In Chapter 3, a new 
technique is developed that uses a reduced-order model of a known larger 
system and spectral assignment procedures to reassign selected eigen- 
values of the system. This is accomplished without affecting the eigen- 
values and eigenvectors not included in the reduced-order model. Sec- 
ondly, a technique is developed that uses Luenberger' s [2] reduced-order 
observer and spectral assignment procedures to implement a constant full 



4 


state feedback design using dynamic output feedback. The assigned 
eigenvalues and eigenvectors of the original system are retained using 
this technique. In Chapter 4, a design philosophy and then a corre- 
sponding design procedure are developed for the new synthesis techniques 
presented in Chapter 3. A software package is developed to facilitate 
the design of dynamic output feedback control systems using this new 
philosophy and procedure. The package is included as a new mode to an 
existing spectral assignment computer aided design program. The use of 
spectral assignment with reduced-order models and reduced-order observers 
in an actual design problem is demonstrated in Chapter 5. Results are 
compared to those obtained by an alternate design procedure. A program 
listing and an example of a computer aided design session are included 
as appendices. 



CHAPTER 2 


SPECTRAL ASSIGNMENT PROCEDURE 

In this chapter a background of spectral assignment theory is pre- 
sented to support the development in Chapter 3. Definitions of eigen- 
values and eigenvectors are given. Then the effect of eigenvalues and 
eigenvectors on the time response of a system is presented. Lastly, a 
characterization of the freedom available in selecting eigenvectors for 
a given eigenvalue assignment using constant state feedback is pre- 
sented. 


2.1 System Eigenstructure and Time Response 
The eigenvalues of an nth order real matrix A are the zeros of the 
polynomial det [XI-A]. The eigenvalues, X_.eC, form a self-conjugate 
set. That is, for each complex eigenvalue x-j there exists a complex 
conjugate eigenvalue X , + ^ = X .*. For each eigenvalue X., there is a 
right eigenvector, v.eC n , that satisfies the equation. 


for i = 1,..., n. If the eigenvalues of A form a distinct set, then 
the associated eigenvectors are linearly independent [11]. Equation 
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(2.1) is written for all A.. and as 

AV = VA, (2.2) 

where V = [V x ... V ] and a = diag.(x 1 ,...,x n ). Since the columns of 
V are linearly independent, V is invertible. Therefore, 

A = VAV" 1 . (2.3) 

Similarly, for each eigenvalue A. there is a left eigenvector, 
u^eC 1 ] that satisfies the equation 

u. T A = X.u i T (2.4) 

for i = 1,..., n. The left eigenvector equation is written for all A., 
and u. as 

U T A = AU T (2.5) 

where U = [ui ... u n ]. For distinct eigenvalues, the left eigenvectors 
are also linearly independent [12]. Hence, premultiplying equation 
(2.5) by (U 1 )" 1 yields 


Substituting for A 


A = (U T f 1 AU T . 


( 2 . 6 ) 


from equation (2.3) into equation (2.6) yields 
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VAV 1 = (U 1 )" 1 AU T . (2.7) 

Premultiplying by U T and postmul tiplying by V yields 

U T VA = AU T V. (2.8) 

Since A is a diagonal matrix of distinct eigenvalues, equation (2.8) 
can only be satisfied if U^V is a diagonal matrix. For convenience 
the eigenvectors are usually normalized so that U^V = I or 

U T = V _1 . (2.9) 


The effect of eigenvalues and eigenvectors on system time response 
is now presented. Consider the linear time invariant system in Figure 
2.1 represented by the system state equations 


x = Ax + Bu (2.10) 

and 

y = Cx, (2.11) 


where A, B, and C are the plant, input, and output matrices respec- 
tively and x e R n , u e R m , and y e R^. The system time response 
is determined by solving the differential equation (2.10). Let a change 
of coordinates be defined by 
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x = Vz. 


( 2 . 12 ) 


Tlie transformed system is 


z = V _1 AVz + V -1 Bu (2.13) 

y = CVz. (2.14) 

Substituting from equation (2.3) into equation (2.13) yields 

z = Az + V Bu. (2.15) 

The solution of equation (2.15) is given by [13] 

z(t) = e At z 0 + l e A(t " T) U T Bu(t) dr (2.16) 

where z 0 is the initial value of z(t) at t = 0. Substituting from 
equation (2.12) yields the time response 

x(t) = Ve At U T x 0 + V l e A(t " T) U T Bu(T) dT. (2.17) 

The first term of equation (2.17) is called the zero input response and 
the second term is called the zero state response. 

Expanding the zero input response yields 
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x(t) = [vi, .... v n ] 


X : t 

e 1 ... 0 


0 . . . 


X t 
. n 



T“ 


Ul 

» 


• 

• 

T 


_ u n _ 


xo 


X,t, T . A A V , . 

v n e 1 (ui x 0 ) + . . . + v ln e (u n ) x 0 


X it 


X t 


[ V nl e (u! T x 0 ) + . . . + v nn e n (u n T ) x 0 J 


(2.18) 


From equation (2.18) the ith component of the statevector is determined 
to be 


x.(t) 




Expanding equation (2.19) yields 


(2.19) 


xi(t) 

• 


1 

• < 

►— < 

1 

X 1 1 . T . 
e (ui x 0 ) 

I 

- < 
h— * 

3 ■ 

• 

x*(t) 
n v ' 

II 


+ . . . + 

• 

• 

V 

nn_ 


The zero state response is expanded next. 
u(t) be a vector of unit step functions, u 0 . 


V 

e 


(u n T x 0 ) 


( 2 . 20 ) 


Let the input vector 
for computational ease. 



11 


This yields 


x(t) = V / e A ^ t_T Vb u 0 dx 


= (Ve At ) / e" AT dT (U T B u 0 ). 


Since A is a diagonal matrix, the integral term is written as 


f -At ^ 

J e dx 


f e X i T dr 


. . 0 


-At 

/ e dx 


1 (e" Xlt -l) 


. . . 0 


0 ... 


-1 (e- Xnt -l) 


Premultiplying equation (2.22) by the diagonal matrix e c yields 


At r -Ax. 
e J e dx = 


- — (l-e^*) 


. . 0 


0 . . . 


JL (l-e Xnt ) 



The vector K is, defined to be 


K = U B Uq . 


(2.24) 


Substituting equations (2.23) and (2.24) into (2.21) gives 


x(t) =[v 1 . 


• v n^ 


-A (i- e Xlt ) 

x. 


. . . 0 


. X t 

0 . . . -A (l- e n ) 

x 


n -K . X .t 

£ V, ( — I) (l-e ' ) 


i=l 


X . 


(2.25) 


where K. denotes the ith element of K. Expanding equation (2.25) 
yields 


*l(t) 


Vll 

(i) (l-e* 1 1 ) 

v ln 

-K X t 

(— ) (l-e n ) 

• 

• 



+ . . . + 

A 

Xn 


(2.26) 


X .t 

The terms e 1 are called the modes of the system. Equations (2.20) 
and (2.26) show that the eigenvalues of the system determine the rates 
of decay of the modes while the eigenvectors determine the contribution 
of each mode to the various states. Thus, the time response of a system 



can be controlled by proper selection of system eigenvalues atnd eigen- 
vectors. 


2.2 Characterization of Freedom in Eigenvector Assignment 
Given the linear, time invariant controllable system with constant 
state feedback in Figure 2.2, the system state equations are written 

x = (A + BF) x + Bv (2.27) 

y = Cx. (2.28) 

Given that constant state feedback is used, Wonham [14] states that an 
mxn matrix F can be found to assign an arbitrary self conjugate set 
of eigenvalues if the system is controllable. Moore [4] characterizes 
the freedom available to assign eigenvectors for an arbitrary self-con- 
jugate set of eigenvalues. He gives necessary and sufficient conditions 
to find a unique real matrix F that satisfies the eigenvector equa- 
tion 


(A + BF) v. = Vj A. (2.29) 

For i = 1, ..., n when B has full column rank. Associate with each 
eigenvalue X. an nx(n+m) matrix S where 

1 A . 




Figure 2.2. System Model with Constant State Feedback 
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and a compatibly partitioned (n+m)xn matrix 


whose columns constitute a basis for the null space of S . Then the 

A i 

necessary and sufficient conditions to find a unique real matrix F 
that satisfies equation (2.29) are: 

1) Vectors v. eC n are linearly independent, 

* ★ 

2) v.j = Vj whenever X.. = x^, and 

3) v e span (It ). 

i A i 

Thus it is possible to assign an arbitrary selfconjugate set of 

eigenvalues and a set of eigenvectors from within the span of It . The 

A i 

null space N. is determined by the selection of an eigenvalue X.. 

A i 1 

The subspace. It identifies the freedom available to assign eigen- 
A i 

vector . 

2.3 Eigenvector Assignment for Real Eigenvalues 
It is first assumed that X^ eR so that v^ e R n for i = 1,... n. 
Equation (2.29) is rewritten as 

(X.I-A) v.-(BF) v. = 0. 


A i 


M 


X. 

l 


(2.31) 


(2.32) 



16 


Since K, is a basis for the null space of S, , 
A i A i 

that postmul tipi ies 1C gives a resulting vector 

A i 

space of S . Therefore, 

A • 

i 


[X.I-A : B] 


Xl 

M. . 
Xi 


K. = 0. 


Expanding equation (2.33) yields 


[X.I-A] N K, + [B] 


M. . K. 

Xi i 


= 0. 


Since v. e span (N. .), then K. determines where 

1 A 1 I 

subspace v.. exists. Hence, 


v 


i 



K i 


It follows from equations (2.32) and (2.35) that 


Fv. = - M. K.. 
■\ X i 1 

By defining w i as 


w i = - M x. K i» 

equation (2.36) is rewritten in matrix form for all 


then any vector K i 
that 1 ies in the nul 1 

(2.33) 


(2.34) 

in the allowable 

(2.35) 


(2.36) 


i as 


(2.37) 
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F[v lt •••» [ Wl , •••> 


(2.38) 


or 


FV=W. 


(2.39) 


Since the eigenvectors are linearly independent, then 

F = W V- 1 . (2.40) 

2.4 Eigenvector Assignment for Complex Eigenvalues 
It is next assumed that x.eC for i = 1, 2 and X..ER for 
i =3, ..., n. Then the first closed loop right eigenvector must satisfy 
the equation 


[A + BF] (v^£ + j v if^ “ ( V RE + ^ V IM^ ^RE + (2.41) 

where the subscript one is suppressed for simplicity. Equating real and 
imaginary parts yields 


[A + BF] v re = v re X re -v im X m 


(2.42) 


and 
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[A + BF] v IM v IM X RE + v R£ X IM . 


The two equations are written in matrix form as 


OrE 1 ^ : x tm^ : 


IM 


'RE 


-v 


IM 

Fv 


RE 


= 0 


and 


[X re !-A : X TM I : B] 


IM 


'IM 

l-V RE 
-Fv 


IM 


= 0 . 


For the case of complex eigenvalues, the nx(2n+m) matrix is 
as 


S x = t X RE I-AiX IM HB ^ 


and a compatibly partitioned (2n+m)xn matrix *K^ is defined by 


K = 
X 




where the columns of constitute a basis for the null space 
Hence, 


(2.43) 


(2.44) 


(2.45) 


defined 


(2.46) 


(2.47) 


°f s x . 
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= 0 . 


(2.48) 


From equations (2.44), (2.45) and (2.48) it is apparent that the vectors 
in (2.44) and (2.45) are contained in the null space defined by K^. 


Therefore 


1 

LU 

Of 

> 

1 



" V IM 

e SPAN 

p x 

Kd 


M x 


(2.49) 


and 


V IM 


N X 

V RE 

g SPAN 

p x 

~ Fv IM 


M X 


(2.50) 


From equations (2.49) and (2.50) it is apparent that the allowable sub- 


space for v RE and v IM is described by 


RE 

r IM 


g SPAN 


(2.51) 


and 
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V 

V 


RE 

IM 


e SPAN 



Combining equations (2.51) and (2.52) yields 


(2.52) 


RE 

r IM 


N, 


e SPAN 


n SPAN 


(2.53) 


N, 


The following characterization of the freedom available to assign 
complex eigenvectors is not developed here but is proved by Marefat 
[15]. Matrixes a and 3 are defined by 


a 



and 


3 = 



(2.54) 


(2.55) 


K and K„ are defined to be matrices whose columns constitute bases 

a p 

for the null spaces of a and 3 • respectively. Matrix y is defined 
by 


Y 


[K 


Kg] 


T 


(2.56) 
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and is defined to be a matrix whose columns constitute a basis for 
the null space of y. A basis for [SPAN(a)r\SPAN(3)] is [(SPAN(a))' 1 ' + 
(SPAN^)) 1 ] 1 where "+" denotes set direct summation and "l" 
denotes orthogonal complementation. Also, a basis for y 1 is a basis 
for y T . Hence 

e SPAN (K ). (2.57) 



A specific vector within the null space of y is defined by postmulti- 
p lying K by a vector x-p Thus 

= K y x T . (2.58) 



Using equations (2.51) and (2.52), x^ and are defined by 


v 

v 


RE 

IM 



(2.59) 


and 





(2.60) 


The left inverses of 
linearly independent. 



Therefore 


exist since the columns of 


are 
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and 


x cl = 

1 

-L 

v re1 


<< 1 
CL 

' 1 


* < 

» — ! 

t_zL_ 


(2.61) 


x c2 = 

1 

aT* 

i 

-L 

1 

< 

70 

m 

1 


1 

>-* 

i 


v !mJ 


(2.62) 


From equations (2.48) and (2.49) it is apparent that the vectors [-Fv R ^] 
and [-FVj^] lie in the space defined by the columns of M^. Hence 


-Fv re - M x x cl 


(2.63) 


and 


- Fv IM = M X x c2" 


(2.64) 


Since X x and x 2 e C and xi e R for i = 3 , .... n, then 
* 

X 2 = X^ because the eigenvalues form a self conjugate set. Further- 
more, the second condition of spectral assignment requires that v 2 = 

★ 

v^ . Thus the specification of one complex eigenvalue and eigenvector 
contains all the essential information of the complex conjugate pair. 

It is also important to note that if v 2 = v 2 and the pair Vj^ , v 2 are 
linearly independent, then v RE and v^ are also linearly inde- 
pendent. In order to calculate the feedback matrix F, the following 
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definitions are given: 

w, = M x x cl , 

W 2 = -M x x c2 , 

Vi - V RE , 
and 

v 2 = v im* 

Recalling that for the case of real eigenvalues 


(2.65) 

( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 


w i " -M Xi x i (2.69) 

and 

V, - »,, (2.70) 

equation (2.38) is rewritten so that 

F ^- v RE’ v IM’ v 3’ v n-* = ^" M X x cl , ' M X x c2 ,_N X 3 x 3’ '••’" M X n x n^‘ ^ 2 ' 71 ^ 


Substituting equations (2.65) through (2.70) into (2.71) yields 

F[v 1# .... v n ] » [w lt •••» w n l (2.72) 

or 


FV = W. 


(2.73) 
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As in the case for real eigenvalues, 

F = WV- 1 . (2.74) 

This development is easily extended to more than one pair of complex 
conjugate eigenvalues. 


2.5 Use of Eigenvector Freedom 

It is shown in the previous three sections that eigenvectors for 
the selected eigenvalues must reside in an allowable subspace that is 
determined by the plant matrix A, the input matrix B, and the se- 
lected eigenvalues X.. . Normally the eigenvector assignment that is 
most desirable for a given set of eigenvalues is not achievable because 
it does not lie within the allowable eigenvector space. In this case it 
is desirable to select the allowable eigenvector that is closest to the 
desired eigenvector. This is accomplished by projecting the desired 
vector into the allowable space so that the error between the desired 
and the assigned vector is minimized in a least squares sense as illus- 
trated in Figure 2.3. 

The desired vector is projected onto the allowable space by the 
projection operator 



"x <vv 


-1 



(2.75) 


for real eigenvalues and 
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P KY = k y (SVV (2 * 76) 

for complex eigenvalues [15]. Indicating the desired vector by the 
subscript "D" and the assigned vector by the subscript "A", the 
projection is accomplished by the equations 


and 



(2.77) 



(2.78) 


for the real and complex assignments, respectively. 

2.6 Improvement of Initial Assignment by Gradient Search 
The selection of eigenvalues and eigenvectors for a system is nor- 
mally motivated by the desire to shape the time response as discussed in 
Section 2.1. However, once the desired time response is approximated, 
there are often other aspects of the assignment that are unacceptable. 

An example is an assignment which requires extremely high feedback gains 
which are expensive to implement and very sensitive to noise. Another 
example is extreme eigensystem sensitivity to small plant parameter 
variations or modeling errors. The freedom available to select the 
eigenvectors often provides a means to drastically improve these second- 
ary design objectives while only slightly modifying the initial eigen- 
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vector assignment and thus the time response. This improvement is ac- 
complished by modifying the eigenvectors within an area local to the 
original assignment. The vectors are modified in such a manner as to 
reduce the undesirable aspect of the assignment most rapidly. 

A cost function J is defined so that a reduction in the value of 
J corresponds to reduction of the undesirable aspect of an initial 
eigenvector assignment. A gradient matrix is computed in terms of J 
to determine how the eigenvector assignment is most efficiently changed. 
Recalling that the eigenvectors are determined by the equation 

»i*Vi (2 - 79) 

for the case of real eigenvalues and 

V iRE 
v i IM 

for the case of complex eigenvalues, it is apparent that small varia- 
tions in X. will cause correspondingly small variations in the eigen- 
vector assignment. A matrix X is defined as 

X = [Xj, X n ]. (2.81) 

Since this matrix designates which eigenvectors are assigned, it is 
called the designator matrix. A gradient matrix [GR] with elements 
[GR]^ is defined to be 


X i 

X i+1 


(2.80) 



The designator matrix X is then varied according to the rule 


X ij (q+1) = X ij (q) ‘ d[GR] ij (2.83) 

where d denotes the step size of X.. during each iteration. The 

* J 

gradient search is continued until a satisfactory compromise between the 
reduction in the value of the cost function and the modification of the 
time response is achieved. 



CHAPTER 3 


SPECTRAL ASSIGNMENT USING REDUCED-ORDER 
MODELS AND REDUCED-ORDER OBSERVERS 

In this chapter, the use of reduced-order models and reduced-order 
observers in the design of feedback controllers is investigated. A 
reduced-order model of a known system is formulated. It is then used to 
design a constant full state feedback matrix for the original full -order 
system. It is shown that the eigenvalues and eigenvectors reassigned in 
the reduced-order model are reassigned in the full-order system while 
those not included in the reduced-order model remain unchanged. The 
constant state feedback matrix is then implemented by output feedback 
using a Luenberger [2] reduced-order observer. It is shown that the 
eigenvalue and eigenvector assignments in the full-order system remain 
unchanged when a reduced-order observer is used. 

3.1 Motivation for Using Reduced-Order Models 
and Reduced-Order Observers 

Models representing the behavior of physical systems often consist 
of a very large nunber of coupled linear differential equations. Such 
models are difficult to use when designing control systems due to exces- 
sive requirements for computer time and memory, and to the numerical 
analysis problems inherently present when dealing with large systems of 
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equations. It is, therefore, desirable to develop a design procedure 
which utilizes reduced-order system models. 

The spectral assignment synthesis methods described in Chapter 2 
use full state feedback. However, full state feedback is not feasible 
for most systems because there are often fewer outputs than system 
states. Full state feedback is implemented by the use of a dynamic ob- 
server for these systems. The use of a full system observer is unneces- 
sary since some states usually are available for measurement and there- 
fore need not be estimated. A reduced-order observer is therefore de- 
sirable in order to minimize the complexity of the control system. 

This chapter develops a reduced-order model and reduced-order ob- 
server. Control system design for the full-order system is accomplished 
using the reduced-order observer. Reduced-order models and observers 
have been used for several years. However, it is shown here that the 
eigenvalues and eigenvectors assigned using the reduced-order model are 
retained in the closed loop full-order system while the eigenvalues and 
eigenvectors not included in the reduced-order model remain unchanged in 
the closed loop full-order system. 

3.2 Reduced-Order Model Formulation 
A reduced-order model that is used in the design of a constant 
state feedback controller for the full-order system model is formulated 
in this section. The reduced-order model contains the eigenvalues that 
are to be reassigned in the full-order model. Let the original system 
model be described by the state equations 
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x = Ax + Bu (3.1) 

and 

y = Cx (3.2) 


where x eR n , ueR m , and yeR* 5 . It is assumed that the eigenvalues 
of A are distinct and are. denoted by (X 1# ...,x n ). The corresponding 
modal matrix for A is denoted by V = [v^, ...» v ] where v. de- 
notes the eigenvector corresponding to A.. The system model is trans- 
formed by defining a new state variable 


-l 

z = V x. 


(3.3) 


Equation (3.1) is transformed to give 


z = Az + Bu (3.4) 

. - 1 , -l 

where A = V AV = diag (\ l , . . . , x p ) , and B = V B. The system is 

now partitioned to separate the eigenvalues to be reassigned in the 

reduced-order model from those that will remain unchanged in the full- 

order system model. Thus, 


» 


l — 
o 
• • 

<r 


V 


"B ' 

Zi 

= 

i 



+ 

1 

• 

Z 2 


0 • A 2 


Z 2 


B 

* 


# 



2 


u 


(3.5) 
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where z^R 1 * and z 2 eR n ” k . The eigenvalues included in the reduced- 
order model must be contained in A x and those which are not included 
must be contained in A 2 . The reduced-order model is thus described by 
the state equation 


z i = Ai zi + Bi U. ( 3 . 6 ) 

The reduced-order model is then used in conjunction with the spectral 
assignment procedure to assign eigenvalues and partially assign eigen- 
vectors in the full-order system model. However, the relationship 
between the eigenvalue and eigenvector assignments in the reduced-order 
and full-order models must be investigated first. 

3.3 Spectral Assignment Using Reduced-Order Models 
The reduced-order model is used to design a constant state feedback 
matrix for the full-order system. The relationship between the eigen- 
values and eigenvectors of the closed-loop reduced-order model and the 
closed loop full-order system must be understood in order to accomplish 
this. The relationship between reduced-order and full -order system 
eigenvalues is determined first. Let F denote a constant state feed- 
back matrix computed for the reduced-order model. The control law is 
then written as 


u = F zi. 


(3.7) 


The reduced-order model closed loop equation is therefore 
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Zi = (Ai + Bi F) Zi. (3.8) 

F is now implemented about the full -order system by assuming full state 
availability in the full order model and transforming F back to the 
original coordinate system. Equation (3.7) is rewritten as 

u = [F:0] = F'z. (3.9) 

L 1 

Substituting for z from equation (3.3) yields 

u = F •V~ 1 x = 7x. (3.10) 

Hence, the closed-loop full order system is written as 


x = Ax + B Fx 

= [A + B F]x. (3 - U) 

The eigenvalues of the full order system are the eigenvalues of 
[A + BF], This matrix is rewritten as 


[A + BF] = [VAV -1 + VBF'V -1 ] 
= V[A + BF'] V" 1 . 


(3.12) 
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Since [A+BF] and [A + BF 1 ] are related by a simil arity transform- 
ation, they have the same eigenvalues. Matrix [A + BF'] is expanded 
as 


[A + BF'] = 

1 

l> 

► |.* • 
lO 

1 

+ 

liEIQ 

A A 


L° 5A2 _ 


B 2 F:0 

f 

m 


_ 


- 


A . 

BlEiO 


• 


A A • 

B 2 F JA 2 


(3.13) 


The eigenvalues of this matrix are obviously the eigenvalues of 
Ca x + Bj F] and [a 2 ] - Thus it is possible to reassign the k eigen- 
values included in the reduced-order model without modifying the (n-k) 
original system eigenvalues which were not included in the reduced-order 
model . 

The relationship between eigenvectors of the reduced-order and full 
order system models is determined next. The eigenvector equation for 
[A +BF'] is written as 


[A.I-A-BF 1 ] v. = 


A.j I-Ai — B 1 F ; 0 

A A , 

-B 2 F :a.-a 2 J 


v ii 

A 

V 2i 


= 0 (3.14) 


which yields the two equations 


and 


[A-j I-Ai -BiF] v^ - 0 


[-b 2 f] v li + [A.I-A 2 ] v 2i = 0. 


(3.15) 

(3.16) 
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If A. is an eigenvalue of [A x + % l F], but not an eigenvalue of 
[A 2 ], then from equation (3.15) it follows that v^ is an eigenvector 

A A 

of + BjF], Since \ . is not an eigenvalue of [a 2 ], [x . I- a 2 ] is 
nonsingular. Thus, from equation (3.16) 

V 2-j - [A. I-A 2 3 B 2 Fv^ . (3*17) 

Therefore, v i is written as 


[_ [ ^- -j I ~A 2 3 B2F. 

Equation (3.18) illustrates that the first k elements of v^ can be 
assigned using the reduced-order model while the remaining (k-p) ele- 
ments are linear combinations of v^. . 

On the other hand, if A^ is an eigenvalue of A 2 and not an 
eigenvalue of [A x + E^F], the matrix [a.I-A^BjF] is nonsingular. 
Therefore, equation (3.15) is statisfied only if 


v Xi> (3.18) 


v 


li 


= 0 . 


(3.19) 


From equation (3.16) it follows that ^ is an eigenvector of A 2 for 
eigenvalue A... Since A 2 is a diagonal matrix, ^ is written as 



0, k., 0, 


0 ] 


(3.20) 
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where is a nonzero constant. Therefore is written as 


V i = V, 


(3.21) 


Let v i be the eigenvector of [A+ BF] corresponding to eigenvalue 
A... The eigenvectors v^ and v.. are related by the transformation 


Vj = V v ( 


(3.22) 


Expanding equation (3.22) for eigenvalues of [Ai + BiF] yields 


v n •. v i2 


V 21 : V 22 t x i I ” A 2^ B 2 F 


+ V^2 ~ 

^21 ^22 ^ I ”"^2 ^ B^F li 


(3.23) 


Often when the reduced-order model contains only the dominant modes 

-l 

of the system, V 1X = V n + V 12 [A^-A;)] B 2 F. Then the top k compo- 
nents of the first k eigenvectors are assigned by choosing v^ as 


5 u =v n V 


(3.24) 
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When the above approximation does not apply, an initial assignment of 
v li is made using equation (3.24), F is calculated, and the error 
between the desired top k components of the first k eigenvectors and 
the actual assignment is calculated. A gradient search procedure is 
then used to reduce this error. 

For which are eigenvalues of A 2 , equation (3.22) expands to 



Substituting for ^ from equation (3.20) yields 



(3.25) 


(3.26) 


Therefore, the last (n-k) eigenvectors in the closed loop full order 
system are the original open-loop eigenvectors. Thus eigenvectors cor- 
responding to the eigenvalues retained in the full order system model 
remain the same in the final closed loop system. 

3.4 Reduced-Order Observer Formulation 
In order to implement the feedback matrix F calculated in equa- 
tion (3.10), full state availability is required. However, if the 
number of outputs p is less than the number of states n, then (n-p) 
states must be estimated. This section parallels Luenberger' s [2] de- 
velopment of a reduced-order observer system to allow the implementa- 
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tion of a full state feedback matrix in such a system. The following 
section shows that not only are the eigenvalues of the observer system 
retained in the closed loop system, but that the original eigenvector 
assignment is not affected by use of the observer. 

The open loop full-order system model described by equations (3.1) 
and (3.2) is 


x = Ax + Bu (3.26) 

and 

y = Cx. (3.27) 

If the feedback matrix 7 is implemented about the system model, then 

u = Fx. (3.28) 


It is assumed without any loss in generality that the first p columns 
of C are linearly independent. A transformation matrix M is defined 
to be 


M = 


Ci 


0 



(3.29) 


where I p p is the (n-p) th order identity matrix. A new state variable 
is defined by 
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z = Mx. 


(3.30) 


Equations (3.26), (3.27), and (3.28) are transformed and written as 


z = A z + Bu, 


(3.31) 



(3.32) 


and 


u = F z (3.33) 

where 


A = MAM" 1 , 

(3.34) 

B = MB, 

(3.35) 

CM"' = [I p ; 0], 

(3.36) 


and 
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F 



(3.37) 


Equations (3.31) and (3.32) are expanded as 


r n 


r~ 

1 


“ - 


m mm 

Zl 


^11 

A i2 


*1 


Bi 


— 





+ 


_ Z 2_ 



a 22 _ 


_ z 2. 


Jz .. 


(3.38) 


y ■ [Ci C 2 ] 




(3.39) 


Therefore the p components of the output vector are the first p 
states of the transformed system denoted by Zi . The reduced-order 
observer must then estimate the remaining (n-p) states denoted by Z 2 . 
Equation (3.38) is expanded to be 


zi = Ah z i + Ai 2 z 2 + 8i u (3.40) 

and 

A/ A/ fSJ /V *-S J A/ 

z 2 - A 2 i Zi + A 22 z 2 + B 2 u. (3.41) 


rsj 

Since z x is available as the output vector y, it can be different- 
iated to generate Zi . Hence, equation (3.40) is solved for Ai 2 z 2 
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which is used as an input to a reduced-order observer to approximate 
z 2 . The proposed observer is shown in Figure 3.1. It is desired that 

0 = z 2 (3.42) 

in order to implement the feedback matrix F. The observer state equa^ 
tion is written as 


0 = E8 + Lz x - LA U Zi + [n-LBju. (3.43) 

The need to differentiate zi is avoided by redrawing the observer 
as shown in Figure 3.2. If w is defined by 


w = 0-LZi , 


(3.44) 


then 


w = Lz - 0 


= 0 . 


(3.45) 


Using equation (3.42) to substitute for 6 in equation (3.45) yields 


w + L ii - z 2 = 0. 


(3.46) 


Calculating each term of the above equation results in the three equa- 
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tions, 


w = E [w + LzJ + [n-LBi]u - LAnZi, (3.47) 


~ IV 


Lzi - LA11 Z\ + LA12 T- 2 . + LBi u. 


(3.48) 


and 


A/ IV 


-Z - -A21 Zi - A22 Z2 - B2U. 


(3.49) 


Substituting equations (3.47), (3.48), and (3.49) into (3.46) gives 


Ew + [EL-A 21 ] Zi + [LAi 2 ~A 22 ] Z 2 + [fi-B 2 ] u = 0. (3.50) 


Using equations (3.42) and (3.44) to substitute for w yields 


[A 2 i] z x + [E-A 22 + LA 12 ] z 2 + [n-B 2 ] u = 0. (3.51) 

Since the input and state vectors are not generally zero, the multiply- 
ing matrices must all be equal to zero for the equation to be true. 
Solving for the last two terms gives 


E = A22 - LA12 


(3.52) 
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and 

ft = B 2 . (3.53) 

Matrix A 21 is generally nonzero also. This indicates that the pro- 
posed reduced-order observer is not adequate. Thus the observer is 
modified by adding A 21 Z : to w and grouping terms as shown in Figure 
3.3. To remove the summer located after the integrator, it is noted 
that 


EL z 1 - C (A 22 - LA 12 )L] z x . (3.54) 

Using equation (3.54), the reduced-order observer is drawn as Figure 

r\j 

3.4. The eigenvalues and eigenvectors of E = [A 22 - LA 12 ] are deter- 
mined by proper selection of L since the eigenvalues of 
[A 22 - LA 12 ] are also the eigenvalues of [A 2 2 - A 12 L ]. Chapter 2 
describes a procedure for selecting a proper to achieve a desired 
eigenvalue and eigenvector assignment. A guideline for selecting 
reduced-order observer eigenvalue locations is discussed in Chapter 4. 

The reduced-order observer is now used to implement the feedback 
matrix F. Expanding the control law given by equation (3.33) gives 

= F i z i + E 2 z 2 . 



U = [Fj: F 2 ] 


(3.55) 
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Using equations (3.42) and (3.44), a substitution is made for zz re- 
sulting in 


u - [Fi + F 2 L] zi + [F 2 ]w. (3.56) 

The control law in equation (3.56) is implemented in Figure 3.5. 

3.5 Effect of Reduced-Order Observer on Full System Eigenstructure 
The effect of a reduced-order observer on the eigenvalue and eigen- 
vector assignment in the closed loop full-order system is developed in 
this section. Luenberger [2] has proven that the eigenvalues of the 
original system assignment and the observer assignment remain unchanged 
in the closed loop full-order system. It is shown here that the eigen- 
vector assignment also remains unchanged. The following definitions for 
matrices G, R, and T are given to reduce the algebraic complexity of 
this development. Let 

G = [(A21 - LAn) + (A22 - LAi2)L], ( 3 . 57 ) 

R = [Fi + F 2 L] (3.58) 

and 


T = [-L : I]. 


(3.59) 






Figure 3.5. Control Law Implemented 



-p* 

kO 
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The eigenvalues of the system are determined first. Using equa- 
tions (3.57), (3.58), and (3.59), the total system state equation is 
written as 


z 


w 


A + BRC BF 2 



TBRC + GC 


E + TBF 2 


LW 


(3.60) 


Equation (3.60) is now transformed to an upper triangular form so that 
the system eigenvalues are apparent. The transformation matrix P is 
defined by 



and a new state variable is defined to be 


(3.61) 


v = w- Tz 


(3.62) 


so that 




p - 





_ _ 


■ rv - 

Z 


I 0 


A + BRC 

bf 2 


T 0 


z 

• 

_v _ 


_-T I_ 


/V A/ 

TBRC + GC 

AiA/ 

E + TBF 2 


T I_ 


-V - 


Equation (3.63) is simplified to give 


(3.63) 
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z 


A + B (RC + F T) 


BF 2 



-TA + GC + ET 



v 


(3.64) 


The submatrix [-TA + GC + ET] is equal to zero. This is shown by eval- 
uating each term within the matrix. The individual terms are given by 



-TA = [LAn - A 2 i : LA 12 - A 22 ] » 

(3.65) 

and 

GC = G[I:0] = [A 2 i - LAn + (A 22 - LAi 2 ) L : 0], 

(3.66) 


ET = [-(A 22 - LA 12 )L : -LA 12 + A 22 ]. 

(3.67) 


Hence, 


-TA + GC + ET = 0. 


(3.68) 


The expression (RC + F 2 T) is equivalent to F. This is shown by ex- 
panding (RC + F 2 T) as 


RC + F 2 T = [Fi + F 2 L] [I : 0] + F 2 [-L : I] 


= [F i *: F 2 ] = 


F. 


(3.69) 
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Therefore the substitution of equations (3.68) and (3.69) into (3.63) 
gives 


Z 


v 


A + BF 
0 



(3.70) 


Thus the eigenvalues of the 'system are those assigned to + BF] and 
[E ] . In other words, the use o‘f the reduced-order observer has no 
effect on the original system eigenvalue assignment. 

The eigenvectors of the system are determined next. Transforming 
the state equation using equation (3.30) results in 


x 



0 




L° 

ij L 

0 


A + BT 

1 

cM 

*IJL- 

CO 


V 

0 

E J 


V 



(3.71) 


Let the closed loop system matrix in equation (3.70) be denoted by Ay. 
An eigenvector of Ay corresponding to X is denoted by Vy. Eigen- 
vector Vy is compatibly partitioned so that 
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The eigenvector equation for Ay is 



Equation (3.73) is expaned giving 


(3.72) 


(3.73) 


[ X I ”Aj] Vy^ ~ [BF 2 ] = ^ (3.74) 

and 

[XI - E] v T2 = 0. (3.75) 

It is assumed that the eigenvalues for [Ay] and [E] are distinct 
since their locations are arbitrarily assigned as discussed in Chapter 
2 . 

-1 

Suppose X is an eigenvalue of [Ay] but not [E]. Then [Xl-E] 

-1 

exists. Premultiplying equation (3.75) by [XI-E] yields the result 

v J2 = 0. (3.76) 

Hence equation (3.74) is simplified to 
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[A I -Ay] V-Q ~ 0. 


(3.77) 


Therefore is an eigenvector of [A.j.]. 

Now suppose A is an eigenvalue of [E] but not [A^]. Then 

-l 

[A 1-Ay] exists. Equation (3.75) implies that v^ must be an eigen- 
vector of [E] . Premultiplying equation (3.74) by [AI-A^-] and rear- 
ranging terms gives 

V T1 = [AI-A r r 1 [lff 2 ] v J2 . (3.78) 


Hence for A that are eigenvalues of [A^.], 


where v.^ 

CE], 



L 0 J 

is an eigenvector of [A T ] . For A 


(3.79) 

that are eigenvalues of 


v T2 (3.80) 

where v T2 is an eigenvector of [E]. Thus it is shown that a reduced- 
order model can be used to design a constant state feedback controller 
for a full-order system. The eigenvalues and eigenvectors assigned 
using the reduced-order model are retained in the full system while the 


V T = 


(A I -A t ) B F 2 



eigenvalues and eigenvectors not included in the reduced-order model 
remain unchanged in the full -order system. It is also shown that a 
reduced-order observer can be used to implement a full state feedback 
design without affecting the eigenvalues and eigenvectors of that 


design. 



CHAPTER 4 


DESIGN PROCEDURE 

The methods described in Chapters 2 and 3 are the basis for devel- 
oping a design philosophy, and then a corresponding design procedure, 
for constant state feedback controller design. The design procedure 
presented in this chapter is most useful when a designer is able to 
characterize the desired system in terms of the closed loop eigenvalues 
and eigenvectors as well as the time response. This chapter reviews an 
existing spectral assignment design philosophy. A corresponding design 
procedure and computer aided design package [11] are discussed next. 

Then an extension of the design philosophy is presented followed by a 
corresponding design procedure. This design procedure is included as a 
supplement to the computer aided design package. Lastly, the signifi- 
cant portions of the additional computer aided design software are de- 
scribed in detail. The modified design procedure uses reduced-order 
models and reduced-order observers with spectral assignment methods to 
reassign selected eigenvalues and eigenvectors in the full-order system 
model . 

4.1 Design Philosophy for Full-Order System Models 
The constant state feedback design philosophy for full-order system 
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models is illustrated in Figure 4.1. The objectives faced by a system 
designer are often many and sometimes conflicting in nature. However, 
the location of eigenvalues and eigenvectors, and the system time re- 
sponse are generally the prime consideration. After these objectives 
are satisfactorily achieved, secondary design objectives are considered. 
These secondary objectives include feedback gain reduction, minimization 
of closed-loop system sensitivity to modeling errors or parameter varia- 
tions, and noise suppression. The spectral assignment design procedure 
achieves a satisfactory control design by selecting an appropriate set 
of eigenvalues and approximating a desired set of corresponding eigen- 
vectors. The eigenvalues determine the rates of decay of the various 
system modes while the eigenvectors determine the relative contribution 
of each mode to the different system states and outputs. After a satis- 
factory time response is achieved with an initial eigenvalue and eigen- 
vector assignment, the secondary design objectives are considered. The 
freedom available to select the eigenvectors often provides a means to 
drastically improve these secondary design objectives while only slight- 
ly modifying the initial eigenvector assignment and thus the time re- 
sponse. This improvement is accomplished by modifying the eigenvectors 
within an area local to the original assignment. The direction and 
magnitude of the eigenvector modification is determined by a gradient 
search procedure as discussed in Section 2.6. 

4.2 Design Procedure for Full-Order System Models 
A computer aided design package written by Marefat [11] currently 
exists and is illustrated in Figure 4.2. The software package consists 




Figure 4.1. Eigenvalue/Eigenvector Assignment Design Philosophy 















60 


of several special purpose subroutines that are accessed by the main 
control program. The subroutines may be entered in any order to imple- 
ment specific design objectives according to the design philosophy in 
Figure 4.1. The system description (i.e.. A, B, C) is entered in Mode 
1. An arbitrary set of eigenvalues is assigned in Mode 2, which then 
formulates the allowable subspaces for the eigenvectors. The desired 
eigenvectors are approximated in Mode 3 by projecting them into the 
allowable eigenvector subspaces. Mode 4 allows the system designer to 
observe the time response for various initial conditions and system 
inputs. The initial design is then improved by alternating between 
Modes 4, 5, 6, and 7 until a compromise between primary and secondary 
design objectives is achieved. Modes 5, 6, and 7 modify eigenvector 
components, reduce feedback gain, and reduce system eigen-sensitivity, 
respectively, using gradient search procedures. 

4.3 Design Philosophy for Reduced-Order System Models and Observers 

The design procedure discussed in the preceding section is useful 
for systems where full state feedback is feasible. Another feature of 
the procedure is that it assigns all eigenvalue and eigenvector loca- 
tions. A system designer is often satisfied with several open loop 
eigenvalue and eigenvector locations in a large system. The reassign- 
ment of the remaining eigenvalues and eigenvectors is better accomplish- 
ed using a reduced-order model that contains only those eigenvalues, due 
to reduced requirements for computer time and memory. Also, large sys- 
tems typically have fewer independent outputs than states. A full state 
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feedback design is implemented in this case with an observer system. 

The observer estimates the system states in order to implement the feed- 
back control law. Since some of the states can be obtained from the 
outputs, only the remaining states need to be estimated with an observ- 
er. A reduced-order observer is desirable in this case. It is designed 
using less computer resources than a full-order observer. Also, less 
hardware is required for implementation of a reduced-order observer. 

A design philosophy that uses reduced-order models and reduced- 
order observers is illustrated in Figure 4.3. In order to design a 
control system using spectral assignment with reduced-order models and 
reduced-order observers, a designer must have knowledge of a desired set 
of system eigenvalues and eigenvectors. The original open loop eigen- 
values and eigenvectors are compared with the desired eigenvalues and 
eigenvectors. A decision is made as to which of the eigenvalues and 
associated eigenvectors are satisfactory and which need to be reassign- 
ed. The spectral assignment design procedure is used to assign the 
desired eigenvalues and approximate the desired partial eigenvector as- 
signment using the reduced-order model. Error between the initial 
eigenvector assignment and the desired eigenvector assignment is then 
reduced by a gradient search. Next, the resultant reduced-order feed- 
back matrix is transformed to the full-order system. 

If all of the states are simultaneously available for measurement, 
then the full state feedback matrix is implemented. However, if some 
states are not available, then a reduced-order observer is designed. 

The eigenvalues of the observer are assigned to be slightly more 
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Figure 4.3. Reduced-Order Design Philosophy 
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negative than the dominant eigenvalues of the closed loop system design. 
This is done to ensure that the observer can respond quickly enough to 
follow the states being estimated. Theoretically if the observer eigen- 
values are assigned to be very large negative numbers then the observer 
will provide a better estimate of the states. However, this is not done 
in practice because the observer then acts like a differentiator and is 
very susceptible to noise. 

This philosophy and the synthesis methods described in Chapter 3 
are used to develop an extended design procedure that exactly reassigns 
an arbitrary subset of the original system eigenvalues which are includ- 
ed in a reduced-order model. A partial eigenvector assignment is then 
approximated for these eigenvalues. This control is implemented with a 
reduced-order observer if there are fewer system outputs than states. A 
contribution of this thesis is that the reduced-order design and imple- 
mentation are accomplished with the knowledge that the eigenvalues and 
eigenvectors not included in the reduced-order model remain unchang- 
ed. 


4.4 Design Procedure for Reduced-Order System Models and Observers 
The computer aided design package written by Marefat has been modi- 
fied as illustrated in Figure 4.4. An additional mode (Mode 9) has been 
added to incorporate the use of reduced-order models and reduced-order 
observers in system design. The full-order system description is enter- 
ed in Mode 1. If a reduced-order model is to be used in the control 
system design. Mode 9 is entered. Otherwise the design procedure con- 
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tinues as described in Section 4.2. 

A flowchart for Mode 9 is shown in Figure 4.5. Several of the 
existing program subroutines including Modes 2, 3, 4, and 6 are called 
from Mode 9. Three new subroutines are also called from within Mode 9. 
These subroutines are described in Sections 4.5 and 4.6. 

After Mode 9 is entered the reduced-order model is formulated. 

Mode 2 is then called automatically and the reduced-order model eigen- 
values are assigned. The designer is now prompted to enter the desired 
partial eigenvector assignment for the full-order system. An initial 
eigenvector assignment is calculated from the reduced-order model using 
equation (3.24) and the assignment is approximated using Mode 3. A 
gradient search is then initiated in order to decrease the error between 
the desired and actual partial eigenvector assignment. Upon completion 
of the gradient search, the reduced-order model feedback matrix is cal- 
culated and transformed to the full-order system coordinates. A re- 
duced-order observer is formulated next. Eigenvalues and eigenvectors 
are assigned to the observer using Modes 2 and 3. Finally a time re- 
sponse is calculated and displayed for the combined system using Mode 4. 
If the designer is not satisfied with the time response Mode 9 is re- 
entered. 

Two portions of the above design procedure required an extensive 
programming effort. Calculation of the cost function used in the gradi- 
ent search is described in Section 4.5 and the gradient matrix calcu- 
lation is described in Section 4.6. 




Figure 4.5. Mode 9 
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4.5 Cost Function 

The cost function J that is used in the gradient search routine 
is a measure of the error between the actual and desired partial eigen- 
vector assignments in the full system model. Calculation of the cost 
function is accomplished in two parts. The actual partial eigenvector 
assignment is computed first using a subroutine called VACT. The actual 
partial eigenvectors are then used to compute the value of J in a sub- 
routine called ROCOST. 

A flowchart illustrating VACT is given in Figure 4.6. The actual 
partial eigenvector assignment is denoted by v^ and the eigenvector 
assigned in the reduced-order model is denoted by v^,.. This is consis- 
tent with the notation used in Chapter 3. The partial eigenvector as- 
signment v^ of the full order closed loop system that is obtained by 
assigning v^ in the reduced-order model can only be determined after 
all reduced-order model eigenvalues and eigenvectors are assigned and 
the feedback matrix F is computed. The top half of equation (3.23) is 
given by 

v ii = t v n + v i2 £V' a 2^ Vhr ( 4 * 1 ) 

For a real eigenvalue the subroutine computes v^ using equation 
(4.1). If the eigenvalue A., is complex the calculation becomes 
slightly more involved. Separating equation (4.1) into real and imagin- 
ary parts yields 
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v li < v ll +v l2( x i I - A 2 ) ’ 1 V^li ' [V 12 (X i 
11 RE 11 'RE * * 11 RE U 1 IM ^ 11 IM 

and 

v li - ^ll + ^12^ X i ^ -A 2^ + ^12^ X i ^ ^V^li 

IM 11 ^ RE * ^ il IM Ld n IM 1 in RE 


m] >N /\« 

-*■ n n 


(4.2) 

(4.3) 


Equations (4.2) and (4.3) are used to compute partial eigenvector as- 
signments for complex eigenvalues. The actual eigenvector assignments 
are then used in subroutine ROCOST to compute the value of J. 

A flowchart illustrating ROCOST is given in Figure 4.7. If the 
desired partial eigenvector assignment is denoted by v^ and the actual 
partial eigenvector assignment is denoted by v, then the cost function 
is calculated by 


J + E < v ii- v n >%i (4.4) 

ij 1J u ij 1J 

where a., are arbitrary weighting constants. The weighting constants 
determine the relative penalty between the eigenvector component errors. 
For example, if one eigenvector component has a much larger weighting 
constant than the others, then an error in that component receives a 
much greater penalty than other component errors. 

4.6 Cost Function Gradient 

The cost function gradient matrix is computed in subroutine ROGRAD. 
A flowchart illustrating ROGRAD is given in Figure 4.8. It is seen from 
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equation (4.4) that the cost functon J is a function of the partial 
eigenvector assignment v. Hence, it is also a function of the desig- 
nator matrix X which is discussed in Section 2.6. By computing a 
gradient of the cost function J with respect to the elements of the 
designator matrix X.^, it can be determined how to vary the designator 
matrix in order to reduce the cost function and therefore the error 
between v Q and v. Recalling equation (2.80), the gradient matrix is 
defined to be 


[GR]^- - 


3J 

ax. . 



|3J ! 

3X ij 


(4.5) 


Solving for 3J/3X. . gives 


3 J 


V pq" V D 

- E 2 ) at Hi) 


«,j P« ” pq D p« «« 


(4.6) 


Since v Q is a constant valued matrix, 


. . 3 V 

— - 2 2 cc (v„ -V n ) pq 


« M M °Pq axT ‘ 


(4.7) 
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To evaluate aVp^/SX^., q is substituted for i and the subscript 1 
is dropped in equation (4.1) to yield 


v q ■ t v n + ^(y-vViv 


(4.8) 


“ ™ t h ~ u _ 

Since v is the p element of vector v , then 3v /3X.. is the 

th ^ q qq ^ 

p element of 37 /3X... Solving for 37/3X.., then gives 

H 1 J M ' J 




SX t 


3x 1j 


3 V 


= V 11^- + [V 12 (x l-A 2 )- l8 2 ][F 


3X ij 


as 

— ?.]. (4.9) 

3 X . . 3X . . q 
1J 1J 


In order to evaluate 3F/3X. equation (2.73) is modified to be 

* J 


F = WV 


-1 


(4.10) 


and the element of F in the p^ row and q^ column is denoted by 

f so that 
pq 


f 

pq 


- R pM[v _1 ] C q . 


(4.11) 
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X. L_ 

Rp is a row vector with a one in the p tn column and zero elsewhere and 

th 

C q is a column vector with a one in the q n row and zero elsewhere, 
so that 


3f 


pq . 


9X ij 


- Rn [ 


3W w_] 


+ W 


9X ij 


9V- 1 

9X ij 


i c«- 


(4.12) 


Solving for 9W/3X^. yields 


9 W 



9X u 


K, 



(4.13) 


From equation (2.68) , 




(4.14) 


Therefore 


au 9X i 

— " [-CM, ] L [-M ] -]. 

9 X • • 1 9X,. a k 

• J * J * J 


(4.15) 


j. ^ 

Noting that only the j column of X is dependent on X^, 

' J 


9U 

— =- = [0, .... 0, [-M. ] 3X./9X.., 0, ..., 
3 X i j 1 J 1J 


0 ]. 


(4.16) 
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Furthermore, only the ith row of X. is dependent on X... It is 

J 1 J 

easily seen that 8Xj/3X.j is a column vector with a one in the i^ row 
and zero elsewhere. Hence, equation (4.16) is written 


3W 

9X ij 


“ •••> o, 0, ..., 0] 


(4.17) 


J. L 

where [-M xi 3 i denotes the 1 column of C -M, - ] . It is noted that 




8V- 


3X . . 


= v- 1 


8 V 
3X. . 


V-l 


and 


3 V _ J 

8 X i j 8X ij 


[vi, ..., vj. 


(4.18) 


(4.19) 


t h 

Since only the j eigenvector is a funtion of X.., it follows that 

* J 


3 V 


3V. 


[o, •••» 0, — — , 0, ..., 0]. 

3 X . . 3X . . 

1J 1J 


(4.20) 


Substituting from equation (2.77) gives 


3 V 

9X ij 


= [0, ...» 0, [N. ] 3X./3X.., 0, ..., 0] 


= [0, ..., 0, [N. ]., 0, ..., 0], 

A j 1 


(4.21) 
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J. L. 

where [N, L denotes the i tn column of [N, ]. To evaluate av_/3X^.. 

Aj 1 A ^ M 'J 

equation (4.21) is postmulti plied by C^. Similarly, to evaluate 
3 v /3 X . , equation (4.9) is premultiplied by R . 

pq i p 

Hence, it is shown that the partial derivatives are computed by 
selecting appropriate rows and columns from the [N^] and [M^] matri 
ces. This reduces the calculation of the gradient matrix [GR] to a 
bookkeeping operation easily implemented in a computer program. It is 
not necessary to numerically approximate a derivative quantity. Sub- 
routine ROGRAD computes the cost function gradient matrix using this 
procedure. 



CHAPTER 5 


DESIGN EXAMPLE 


The design procedure described in Chapter 4 is illustrated in this 
chapter by an actual design problem. A controller is designed for the 
lateral axis model of an L— 1011 aircraft using a reduced-order model and 
a reduced-order observer. The resulting design is then compared to an 
output feedback controller designed by Andrey et al. [15]. It is shown 
that the design procedure presented in this thesis is a viable tool for 
constant feedback controller design. 

5.1 Original Lateral Axis Model 

The lateral axis model of an L— 101 1 aircraft is used as the orig- 
inal full-order system model. The state vector x is given by: 

xi = r = Yaw rate (Radians/second) 
x 2 = 6 = Sideslip angle (radians) 

X 3 = p = Roll rate (radians/second) 
x 4 = <j, = Bank angle (radians) 

x 5 = 5 r = Rudder deflection (radians) 

Xe = = Aileron deflection (radians) 

x 7 = f = Washout filter state. 

W 
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Rudder and aileron deflections (states 5 and 5) .produce changes in the 
yaw rate, sideslip angle, roll rate, and bank an£)le (states 1-4). The 
coordinate system is illustrated in Figure 5.1. * Under certain con- 
ditions yaw rate is equal to the derivative of the sideslip angle with 
respect to time while roll rate is equal to the derivative of the bank 
angle with respect to time. The washout filter ,is a high pass filter 
for the yaw rate. 

The A, B, and C system matrices are given by 


-0.154 

1.54 

-0.0042 

0 

-0.744 

-0.032 

0 

-0.996 

-0.117 

-0.000295 

0.0386 

0.02 

0 

0 

0.249 

-5.2 

-1.0 

0 

0.337 

-1.12 

0 

0 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

-20.0 

' 0 

0 

0 

0 

0 

0 

0 

-25.0 

0 

0.5 

0 

0 

0 

0 

0 

-0.5 


B 


0 

0 

0 

0 

20 

0 

0 


0 

0 

0 

0 

0 

25 

0 


ri 

0 

0 

0 

0 

0 

-11 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

.0 

0 

0 

1 

0 

0 

0. 


The system input, u, consists of components 




Figure 5.1. Aircraft Axis System 
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and 


u x = 6 rc = Rudder command (radians) 

u 2 = 6 = Aileron command (radians). 

^ dC 


The open loop eigenvalues of this system are: 


A i ,2 = -0.08819 - j 1.269 

X 3 ' = -1.085 

A 4 = -0.00965 

A 5 = -20.0 

A 6 = -25.0 

Ay = -0.5 


- Dutch roll mode 

- Roll subsistence mode 

- Spiral mode 

- Rudder mode 

- Aileron mode 

- Washout filter mode. 


The open loop system time response is shown in Figures 5. 2-5 . 8 for zero 
input and an initial condition of <j>(0) = 1 degree. After ten seconds 
the system states are still oscillating and the bank angle <|> has not 
yet reached zero degrees. 

It is known that a desirable eigenvalue assignment for the system 
is 


and 


A j ,2 = —1.5 - j 1.5 
A 3 = -2.0 ± j 1.0. 


When the roll subsistence mode A 3 and the spiral mode A 4 are a com- 
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Figure 5.8. Washout Filter-Open Loop Response for $(0) = 1 
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plex conjugate pair they are collectively referred to as the roll mode. 
It is also known to be desirable for the roll and dutch roll modes to be 
decoupled. This decoupling is accomplished by the eigenvector selec- 
tion: 


1 


X 


0 


0 

X 


1 


0 


0 

0 


0 


1 


X 

0 

X 

, V 2 = 

0 

X 

• V 

X 

X 

* V4 = 

1 

X 

X 


X 


X 


X 

X 


X 


X 


X 


where X denotes "don't care."' Andry, Shapiro and Chung [15] closely 
approximate the above eigenvalue and eigenvector assignment for this 
system using constant output feedback. Eigenvalue/eigenvector assign- 
ment techniques are used to design the constant output feedback matrix 

3.35 -0.159 -4.88 -0.379 

K 1.42 2.38 -6.36 3.8 

The closed loop time response is shown in Figures 5.9-5.15. The closed 
loop eigenvalues of the design are 

X 1>2 = -1.052 ± j 1.497 

X 3>4 = -2.001 ± j 0.9995 

X 5 = -17.05 

X 6 = -22.01 

X7 


-0.6989. 



89 


TIME 


SI MW* 


R 

E 

S 

P 

0 

N 

S 

E 


* * 

* * 



* OOE+OO .40E+0X •80E+01 *X2E+02 


f I M E 

Figure 5.9. Yaw Rate-First Closed Loop Response for <j>(0) = 1° 





n to 2: o “3 to m 


90 


time si mu* 


* 
* 
* 
* 
+ 
* 
* 
* 
* 
+ 


+ #* 


•60E-02 


-02 


* 

# 

# 

* 

+ 

* 

* 

% 

*. 

+ 



•00E+00 


« 40E+0I 


♦ SO Ilf Oi 


♦ 12irf0 


TIM E 

Figure 5.10. Sideslip Angle-First Closed Loop Response for <j>(0) = 1° 


*«•****-}-*****+*****.-+-*** #■-+-•?**<• ie- *-*-.**•*** 



n o z □ "? w n 


91 




* 

X 

X 

X 



T X H EL 

Figure 5.11. Roll Rate-First Closed Loop Response for <j>(0) = 1° 



U CO C- O (O UJ 


92 


♦20E+01 


« jhe+oi 


•16E+01 


« 14E+01 


♦ 12FL+03 


•30E+01 


+##**** 

* 

* 

* 

* 

+ 

* 

* 

* 

* 

+ 

# 

* 

* 

* 

+ 

4: 

* 

* 

* 

+ 

* 

* 

* 

* 


TIME SJMIJ* 


* 

* 

# 


*X 

*a 

* IX 

•80E+00 + X 

* X 


«60Ef00 f 
* 


♦ 40E+00 + X 

* X 

* k 

* x 

* x 

♦ 20E+00 + X 

* \ 

* XX 

* XX 

% XX 

• OOE+OO +:m*:ms*+fc xXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX $ft*3:t;i 
♦ OOE4 00 •40E+01 «80E+0& 


T I M E 


Figure 5.12. Bank Angle-First Closed Loop Response for $(0) = l e 
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The first four components of the first four eigenvectors are: 


V 


1,2 


1 

0.03066 ± j 0.3488 
-0.0036 ± j 0.0004 
0.0013 ± j 0.0011 



-0.0029 * j 0.0012 
0.0045 ± j 0.0053 
1 

-0.3999 + j 0.2000 


It is noted that by using constant output feedback that the four eigen- 
values X x - x 4 are placed almost exactly and that the roll and dutch 

roll modes are decoupled. However, the other eigenvalues (X 5 - X 7 ) are 
also moved by the design. The design procedure described in Chapter 4 

is now used to formulate an alternate design that exactly places Xi - x 4 
without changing x 5 - x 7 . Eigenvectors for \ l - x 4 are also assigned 
to achieve roll and dutch roll mode decoupling. This is achieved with- 
out modifying the eigenvectors associated with X 5 - x 7 . Furthermore, 
this design is done using a reduced-order model to specify a constant 
feedback matrix for the original full-order system. The full state 
feedback matrix is then implemented by dynamic output feedback. 


5.2 Reduced-Order Model Design 

Since the open loop values of x 5 - x 7 are known to be desirable, 
and the rudder, aileron, and washout filter states are unspecified, no 
reassignment of these modes will be made. Therefore, they are not in- 
cluded in the reduced-order model. On the other hand, \ x - x 4 are to 
be reassigned and are included in the reduced-order model. The full 
order system matrices are transformed by equation (3.3) and partitioned 
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as in equation (3.5) to yield the reduced-order model system matrices 


- 





- 

1 

-0.08819 

1.269 

0 

0 


-1.706 

-0.02580 

-1.269 

0 

-0.08819 

0 

0 

-1.085 

0 

0 

II 

03 

0.3961 

-0.2772 

0.06988 

-0.2878 

0 

0 

0 

-0.009165 


-0.2698 

-0.1528 

- 



- 


_ 



Spectral assignment synthesis methods are then used to assign the eigen- 
values 


- -1.5 - j 1.5 
and 

^3 = —2.0 ± j 1.0. 


A partial eigenvector assignment that achieves roll and dutch roll mode 
decoupling is given by 


“ 


- 

20 


0 

6 ± j 7 

V - 

0 

0 

* V 3,4 

20 

0 

1+ 

00 

1 

. 


L J 


An initial attempt is made to assign the eigenvectors using equation 
(3.24). The partial assignment in the full order system is found to be 


“ 


- ” 

19.44 ± j 0.33 


-0.10 + j 0.04 

6.76 ± j 7.25 

W 

0.08 ± j 0.10 

1.06 + j 0.10 

’ V 3,4 

20.20 ± j 0.11 

-0.42 ± j 0.64 

-8.54 + j 3.91 
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The gradient search routine described in Section 2.6 is now used to 
improve the initial vector assignment. Elements of a weighting matrix 
are entered into the computer and a value is calculated for the cost 
fucntion J using equation (4.4). A cost function gradient is calcu- 
lated as in Section 4.6 and the initial eigenvector assignment is varied 
to reduce the cost function. The weighting matrix is varied to increase 
or decrease the relative importance of each eigenvector component and 
the gradient search is continued. This procedure is repeated until a 
satisfactory improvement of the initial assignment is achieved. In this 
example the final partial eigenvector assignment in the full system 
model is given by 



19.45 ± j 0.34 
6.76 ± j 7.25 
0.45 ± j 0.33 
-0.07 ± j 0.68 



-0.10 + j 0.04 
0.08 ± j 0.10 
20.20 ± j 0.11 ‘ 

-8.54 + j 3.91 


The vectors are scaled to give 


1 


-0.005 + j 0.002 

0.35 ± j 0.37 
0.02 ± j 0.02 
0.003 ± j 0.04 

• V 3,4= 

0.004 ± j 0.005 
1 

-0.424 + j 0.191 


It is seen from the above vectors that the roll and dutch roll modes 
have been significantly decoupled. The required gain matrix in the 
reduced-order model is given by 
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F = 


1.319 

-3.854 


-1.650 0.169 

0.583 -5.810 


-1.724 
31.18 J* 


The constant state feedback matrix in the full order system is computed 
using equations (3.9) and (3.10) to be 

-r _ [3.66 -3.13 -0.176 -0.372 -0.137 0.003 o] 

~ [1.54 -6.03 2.70 4.35 -0.008 -0.120 Oj * 

In order to implement this full state feedback matrix, an observer is 
now designed by the procedure described in Section 3.4. The observer 
eigenvalues, X Qi , are selected so that 


X 01 ” " 5 ’ X 02 " " 6 ’ X 03 " 7> 

This selection makes the observer modes faster than the modes contained 
in the reduced-order model. The observer eigenvectors, Vq^ , are arbi 
trarily assigned to be 



1 


0 


0 

< 

0 

II 

0 

0 

’ V 02 " 

1 

0 

* V 03 

0 

1 


The observer matrices are then calculated to be: 


E = 


-5 

0 

0 


0 0 
-6 0 
0 -7 


9 
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■[i 


G = 


and 


1.287 

1.184 

-0.029 

0.710 

-28.3 

-17.1 

-0.372 

4.35 

-89.2 

-42.2 

-5.29 

2.35 

-83.7 

-0.035 

-17.8 

57.6 

48.5 

0.124 

-0.116 

0.270 

TB 

= 

20 

0 

0 

o" 

25 . 
0 



]■ 


The observer is now used to implement the full order system feedback 
matrix F. The closed loop time response is shown in Figures 5.16-5.25. 
It is shown that the response of the yaw rate and sideslip angle for 
this design are more desirable than for the previous design since there 
is less disturbance and faster settling time for both states. On the 
other hand, the roll rate and bank angle responses are almost identical 
for both designs. The controlling surfaces and washout filter states 
are all well within the physical limitations of the system. This illus- 
trates that a viable constant full state feedback control system can be 
designed for a large system using a reduced-order model and that this 
feedback design can be successfully implemented with a reduced-order 
observer. 


5.3 Summary and Concluding Remarks 
A new design procedure for the control of large systems using re- 
duced-order models, reduced-order observers, and spectral assignment 
techniques is presented. A reduced-order model is formulated containing 
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TIME SI MU ♦ 
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Figure 5.16. Yaw Rate-Second Closed Loop Response for <j>(0) = 1° 
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UnE SXMU 4 
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Figure 5.17. Sideslip Angle-Second Closed Loop Response for 4> (0) = 1° 
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T ] MK SI MU ♦ 


* 

* 
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: + £****«« + 
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Figure .5.18. Roll Rate-Second Closed Loop Response for <j> (0) = 1° 
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TIME 8XMU. 


R 

E 

S 

P 

0 

N 

S 

E 
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Figure 5.19. Bank Angle-Second Closed Loop Response for <}>(0) = 1° 
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TXMt: SXMU. 




••40F.-01 


♦ OOE+OO 



■ ♦ 40FL— Oi 


-.80E-01 


-.12F-+00 


- * X6E.+00 


-.20E+00 


-•24E+00 


— ♦ 28E+00 


♦ OOE+OO «40E+0£ *B0E+0:i. .X2E+02 


T X M E 


Figure 5.20. Rudder Deflection-Second Closed Loop Response for <j>(0) = V 
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Figure 5.21. Aileron Deflect ion -Second Closed Loop Response for <J>(0) 
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TIME 8XMU. 
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Figure 5.22. Washout Filter-Second Closed Loop Response for 4> (0) = 1° 
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TIME H X Mi J 4 
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Figure 5.23. Observer State #1 
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TIME SI HU « 
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Figure 5.24. Observer State #2 - Closed Loop Response for <{>(0) = l c 
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Figure 5.25. Observer State #3 - Closed Loop Response for cf> (0) = 1° 
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eigenvalues of interest from an original full-order s'ystem. A constant 
state feedback matrix is designed for the reduced-order model that, when 
implemented about the full-order system, reassigns the eigenvalues con- 
tained in the reduced-order model while those eigenvalues not included 
in the reduced-order model are retained in the full-order system. It is 
then shown that the full state constant feedback matr*ix for the original 
full-order system is implemented by a reduced-order observer if all of 
the system states are not simultaneously available. 
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APPENDICES 


Appendix A contains a software listing of the modifi- 
cations to the spectral assignment computer aided design 
package discussed in Chapter 4. This is followed by an 
example of an interactive design session in Appendix B. 



APPENDIX A: SOFTWARE LISTING 


SUBROUTINE HODKV 

C FUNCTION! MAIM ROUTINE. FOR COMPONENT MOD 11“ IPA T1UH 
IWHSL ROU1 INFS CALLED! (ftfRGEf ,1/Sl/FH, 

C- SPECTRAL ASSIGNMENT ROUTINES I RUGRAD# KOCOST , SEARCH # DSPLAY »VACT • 

o-LOGicAL devices) input ijnxti s imrptn imrn S 

c - STORAGE UNIT (S) I 111-20# XU»20F.t TAPE .1*1 ,NS, (01-20 1 MS 1 1 

INTEGER OR 100# XI.R( 10) 

COMPLEX UE10( 10) ,Z( 10, 1(>) 

REAL LRURO ( 1 0 > , I. IORO ( 10 > , VO ( 1 0 , 1 (i ) , VI’S C 1 0 , 1 0 ) > VDES < 1 0 • 1 0 ) 

REAl Dl AH:f<10#10)#VT( 10)»V0IRV(10# 10) # Al AM! 10# 10) # BI_AH( 10, 10) 
RE At. ALAM 1 1 ( 1 0 » 1 0 > # At. AH 1 2 < 1 0 # 1 0 ) . AI.AH2 1(10,10)# At. AH 22 ( 1 0 # 1 0 ) 
REAL h< 10# 10) tKTILl ( 10# 10) »FTIl.2(10# 10) #1.(10,10) ,K0BS(10, 10) 
REAL DL AH 1(10,10). EL AH1 (10,10)# Cl. AH2 ( 1 0 # 1 0 ) , El. AH (10,10) 
C-RANDOH ACCESS FILES! SYSTEM. FORXX WHERE XX»20L.I , CURRNT 
REAL At. (10# 10) ,11(10, 10) #ELAH( 10,10)# MINV( 10, 10) 

REAL XX(10#10), VA(20),E(20),X(20)#LRE(10) ,I.IH( 10) # UJ( 10) 

REAL U(10,10),V(10#10),VIK'V(10#10)#r(10,10)# AHAT (10,10) 

REAL A( 10, 10) ,»( 10, 10) ,C(10,10) #UKAREA( 130) 

LOMMON/HYfi/A # B # C , ZERO # I Util , NS . N I # NO 

common/ai ib/f # aha t/e j t. ih/i-ar/al/ar/o 

COMMON/ VEf;/ V A # E # X # W J # W # XX # V , V I N V 

COMHON/Rn/OR 100 # VO # VFS # VDES # BLAH? # I.RORO , I. IORO # ALAH22 
DI HEMS IUII CHAR (2) 

EXTERNAL ROCOST , HOUR AO# HOPE? » L10DE3 # MODE 4 
CALL IJEKHETCI.LEVOI.D) 

upewfile « 'roimfo' , access* 'direct * , RECt.o 1 02 # uni r-3?) 

111-20 

READ ( Ilf , REC«1 ) NS # NI , NO, I DOT , ZERO 
REALM X1J#REO2)((A(II,IJ)fI.f'*lfNS)#XI‘*lrN0) 
READ(IU,REC-;i)((B(II,I.I)#I.I-l,NI),XI-l#NS) 
READ(X0#REC«4)((C(II,X.I),X.I»1#NS)#II-1*N0) 
cmmmiMiti enter original eigenvectors ********* *********** 

WRITE(A,H) 

0 FORMAT! IX, 41H WANT TO ENTER NEW ORIGINAL EIGENVECTORS/) 

READ (IS,* IKK 

IF (KK. LE.OXiU TO Al 

J-l 

7 URITE(A# 5) J 

S FORHAT < 2UIIENTER ORIGINAL EIGENVECTOR V,I2) 
READ(t5,*)(V0(I,.l),VT(I),I-l,NS) 

CHTEST«0. 

DU 40 I«1,NS 
CHTESr-CMlEST) VT(I)**2 
40 CONTINUE 

I F ( CHTEHT • l.E • 0 • )00 TO S2 
DO V 1=1# NS 
VOdr.lTD^VTU) 

A CONTINUE 

.(».m 

52 J*JU 

IF ( J*LK«NS)UO TO 7 

WRITE (32# REOI ) ( ( VO ( 1 1 • I J ) , 1.1 s 1 1 NS ) r 1 1 « 1 # NS ) 

Al CONTINUE 

REALM 32# REO 1 ) ( (V()( II, I.)) , 1.1=1 ,N0) , I X' J ’l *NS) 

A1S WRXTE(32.REC-i)((V0(II,IJ),I.^l#NS>.n-l#NS) 

CALL IISVFM (VHV MATRIX! ,V,VO# 10#NS#NS# 4) 

WRITE (A# A 12) 

A! 2 FPRHAl ( 1 X ,31|IWISH TO CIIANDE ANY VALUES OF V?> 

READ(fi,t)KK 
II (KK.LE.OMM TO A13 


UNITE (A, AM) 

414 I- OR HAT < 1 X , 22HENTER I,J, ft NEW VALUE) 

READ(S,*>X,.|,VO(I » J) 

CiCi TO Alt! 

C************ TRANSFORM SYSTEM BY MODAL KATA IX ********************* 
413 CALL LINV2F(V0rNS,l<>#V0XNVf 1DGT, UKARFAr ILK) 

C PAUSE 1 

C CALL USUFM ( Al IVO I N V 1 , A , VO I KV , 1 0 , NS , NS , 4 ) 

call vmji.Kr(vniNv#A,Ns,Ns,Ns,io#io,u, io,^er> 

CALL VMI ILFF (U, VO, NS , NS , NS , 1 0 , 1 0 , Al. AH 1 1 0 # I ER ) 

BO Al A 1*1, NS 
Bli 61 A J*1,NS 
AHSAl.= ADS ( Al. AH ( I r J ) > 

IF < ABSAL . I-T . ZENO > ALAM ( I , J)*F LOAF ( 0 > 

Al A CONTINUE 
C PA OSt 2 

CALL IfSWEH < SHAL AM 1.5, ALAM 1 1 0 « NS t NS , 4 ) 

CALL VMULF'F ( VOINV, B,NSi NSiNI , 10 # 10# BI.AM# 10 , XER ) 

r; pause ;i 

C CALL USUFM ( 15HBLAH 1 1 5 , BLAH . 1 0 , NS , N I , A ) 

CALL VMULFF ( C , VO # NO # NS # NS r 1 0 , 1 0 # CL AH , 1 0 , I Eli ) 

0 PAUSE A 

C CALL 1JSWFM < St 1C LAM IrH. CLAM , 10 , NO , NS , 4 > 

C ************* PARTITION TILDA SYSTEM ****************************** 

BO AOO I«1 , NO 
BO AOO .1*1 , NO 
ALAHU (I,J)»ALAM(I, J> 

AOO CONTINUE 

BO syy 1*1 , NS 
DO 599 »l ir 1 , NS 

I F < ( I . OT . NO > • AND , < .1 . (IT • NO ) > ALAM22 ( I -NO , J-NO ) «r ALAM < I r J ) 

599 CONTINUE 

DO A01 1*1 i NO 
BO A01 ,1*1 ,NI 
BLAMMl t J)irBLAH(Xf J) 

A 01 CONTINUE 

BO A 02 1*1 ,N0 
HO A 02 .1*1, NO 
CLAH1 ( I r J)*CL AH( I , J) 

402 CONTINUE 

C ************* DISPLAY RO MODEL ******** ****************************4 

WRITE (At 004 ) 

004 FORMAT (IX, 25HUISII TO DISPLAY RO MODEL.? ) 

REALMS,* )KK 
IF (KK«LE,(>)00 TO HOS 
UR 1 TE( A, A03) 

A 03 FORMAT! IX,/// ,2011 REDUCED ORDER MODEL,//) 

CALL USUFM( VIIMATRXX Al ,V,ALAHli,10.N0,N()r4> 

CALL USUFM (yUHATRI): Bl ,y,BLAMl,10,N0,NIi4) 

CALL USUFM( VMHATRIX C I , VtCLAMl , 10, NO, NO, 4 ) 

C* *********** TEMPORARILY STORE OR I (i INAL SYSTEM ********************1 
SOS UR I TE ( 32 , RFC* 2 ) NS , N I , NO 

WRITE! 32 rR£C*3) ( (A( XI 1 1,1) , X«I a l,NS) , 11*1 ,NS) 

WRITE ( 32 ,REC»4) < (B( II , X,l> , X J »1 ,NX ) • XX *1 ,NS) 

URITE(32#REC=S) ( <C ( II » I.J) , I,J~ 1 , NS) , XI*i,NO) 
f; ********* STORE RO MODEL INTO FILE *20 ************************** 



Ul< I Tl' ( i!0 * RE1> 1 ) HO » N1 « K1) » I Inn i ZERO 

MRm<iR>tREI>2» < (AUiMtl < . IJ-1 .NCi> . »■*» »NU) 

URnEI20*REC-3M<BLAMiai.IJ>.lJ»»*Mn.I*-t«HI)> 

URnF<2O.REC-4M<i:LAH»(II»IJ*rrJ"lrHO)rII“».R0) 

(IttUI 00 TO H0HE2 AMD AS8IIUI EUif.NVAl.OES EOR RO SYSTEM **»«**»* 

604 FORHAYOX. »»< 111* > • SAII REDUCED ORDER EIOPIVAI.UE ASSI0M1ENT . I7« 1H*> 

call «ni>i:s 

C****«*t ENTER DESIRED PARTIAL EIGENVECTOR ASSIGNMENT ********* 
URITE(6.605) 

605 r ORMAT (1X.44ILENTER DESIRED PARTIAL EIGENVECTOR ASSIGNMENT) 

.1*1 

60V URITE(6»606 >.) 

606 E ORMAT ( IX. 13HEIUEN VECTOR V.I 2) 

READ(t>»*> (VDEG( X » J) » VT( I > » 1*1 .NO) 
if<ads(lim(J>).lt.zl:ro>oo to aoo 
no 607 I-l.NO 

VDESU. JH)*VT(I> 

607 CONTINUE 
J-Jll 
600 J=J11 

C***#*«i* ,, ‘KAfCUEATE f ' AU d’dISPLAY INITIAL 0UES8 *»**«****»»»**»»»*«* 
DO 610 I-l.Nn 
net 610 J-l.NO 

tHi.J>«voa».n 

6 1 ( CALL 1 L I NV2F (O.NO. 10> AHAT . I DOT • UKAREA. I LR ) 

C PAUSE H 

C CALL IJSUFM ( 7HV1 1 INV 1 » 7 . AHAT . 10 # NO . NO. 4 ) 

CAL I. VMUL.FF ( AHAT * VDE8 » NO .NO . NO. 10 .10.0 v 10. XKR ) 

611 rORHAfaxlimiRE THE KOLLOWIHO V MATRIX FOR INITIAL ASSIGNMENT#/. 
12911 RE MEM RE R WHICH V ARE COMPLEX) 

CHARO )«• INITIAL OU* 

CHAR(2)«*E8S FOR VI* 

CALL USWKM < CHAR » 20 # 0 r 1 0 # NO » NO » 4 ) ™™********.* 

C ************ ASSIGN INITIAL EIGENVECTOR OtJESS ********************* 

CALL MODEL! 

XU-20 

REA)) <HI#REI>1) NS.NlfN0.ID0T.2ER0 

READ (XU.RE02) <(A(II.I«I). I*! 31 ! .NS) . II* 1 .NS) 

READ (NJ.RF03) <<B(XI»U>»XJ-itNI)#II-ltN3> 

RKAD(IU.RE(>4>C(G<II.IJ>.IJn.N3>rII-l.N0) 

C PAUSE 6 

C 010 FORMA faX.T/HTME FOLLOUINO IS RO MODEL AFTER MODEJ) 

C CALL IJSUFM < 2HA I . 2 . A . 1 0 . NS . NS . 4 ) 

C CALL USUFM(2HDI.2.Bil0.NS.NI.4) 

C PAUSE 7 

C CALL IJSUFM ( 2HC 1 . 2 . C . 1 0 1 NO . NS » 4 ) 

HO 10 .1-1 . NS 
IU-20W 

OPEN < ACCESS* ' DIRECT' >REGL"202 
1 p UNIT* 10) 

READ ( III.REC«1 ) LRE(J).LIH(J> 

10 CONTINUE 

C PAUSE 0 

DU 012 1-1 .NS 

c WRITE (6.01 l)I.l.m:m.LIM( I ) 

011 FORMAT (IX. 6HL.AMBDA . X 2 . 3X . 5HREAL « . F6 • 3 » 3X • SHI MAO* v F 6 • .1 ) 


012 CONTINUE 

IUT-20INSH 

OPEN (F1LE»m:LIRRNT ' . ACCESS- 'DIRECT' .RECL«iO:2 
1 , UNIT- IUT) 

READ (lllTrREC»l) ((V(XX«1J). X .1*1 .IIS) • 11*1 .NS) 
READ <nif.Ri:C«2> ( (XX( II » I J> . I J : *l . NS) .11=1 »NI ) 
READ (IIJT.REC-4) < <F< II . I J) . X.l*l .NS) . 11*1 .NI ) 
READ ( I LIT . REC«5 ) ( ( AHAT 1 1 1 . 1 J) . 1 J»1 » NS > . 1 1 « 1 . NS ) 
f; CALL IJSUFM ( 10IIMATNIX V I * 10. V. 10. NS. NS. 4) 


C CAL.L IJSUFM < 10HMATKIX XX i r 10 . XX. 10. NI .NS 1 4 ) 
C PAUSE V 

CHAR ( 1 )*M)AIM MATH! * 

CIIAR(2>’**X Ft* 

CALL USWFM (CHAR. 14. Ft 10. Ml .NS. 4) 


C PAUSE 10 

CHAKd )* ‘MATRIX AHA* 


CHAR (2)**TI* 

C CALL IISWFM (CHAR. 12. AIIAT# 10. NS. MS. 4) 

read ( ;i2 . rei>2 > or i so .nix. nox 


X2-ORIOONS 

C********* ASSIGN UEiniUIHO CONSTANTS FOR GRADIENT SEARCH 
VOS DO 70 1*1 .NS 
DO 70 J-l.NS 


******* 


AL ( I . J ) «FL OAT ( 1 ) 

70 CONTINUE 

URllE (6 r VV > 

VV F ORHAT (lX.2£>(lH*)r 27H RO EIGENVECTOR IMPROVEMENT. 20 ( 111* ) ) 

UR I TE ( 6 . V ) 

V FORMAT (IX. 20H WEIGHTING CONSTANTS) 

CAL L USUFI4 < OH WE 1 OHTS I . H . AL . 1 0 . MS , NS . 4 ) 

UR11E(6i 11) 

11 FORMAT ( IX. 1 AH WISH TO CHANGE?) 

READ (5.*> KK 
IF (KK.LE.O) 00 TO 30 
WRITE (6.3) 

S FORMAT (IX. 17HKNTER NEW VALUES I ) 

READ (5.*) ( (AL (I. J).J»l.NS)f I*1.NS) 

C******** CONDUCT GRADIENT SEARCH ******************************* 
30 CALL RGCGGT(CJ) 

WRITE (6.4) CJ 
4 FORMAT (iX.SHC08T«.E1S.6) 


CALL ROGRAD 

CALL BE ARCH (CJ.RGCOST .ROGRAD. S> 

WRITE (IIJT.REOl) •( (VCXIflJ) .XJ^l.RS) »II«1 »NS) 


V04 


C 

902 


C 

C 


WRITET <IilT»REC«2> < ikXdl. IJ) *f?T) 

WRITE (IIIT.RE04) ( (F( II. I J) . XJ»1 .NS) . II-l . NI > 
WRITE (IllT.REl>S> < (AHAT (XI.XJ).IJ ir l »NS> .II-i.NS) 
CALL USWFM CIOIIMATRIX V 1 .10. V. 10. NS. NS. 6) 

CALL VACT 

CALL USWFM (10HHATR IX VKS. 10. VFB. 10 . NS. NS. 4 ) 
WRITE(6.V04) 

FORMAT (IX. 24HUANT TO CONTINUE SEARCH?) 

READtS. t )KK 
If(KK.GT.O)GO TO VOS 


WRITE (6.V02) 

FORMAT (1X.44HU1SH TO DISPLAY 
READ (S.*> KS 
IF (KS.LE.O) GO TO V03 


THE NORMAL. 1 7ED EIGENVECTORS?) 



C CALL DUPLAY (NS .ZERO) 

703 CONTINUE 

C CALL USUFH (lOllflATftXX XX! # 10, XX# l(i»Nl #NS# 4) 

CHAR<1>« 4 (SAIN HA' 

(■HAT«2> = 4 TKIX Ft 4 

CALL imWFH <CHAR,14#F. 10,NX,NS#A> 

CHAIM 1)** 4 MATRIX AHA 4 

CALL El (>RF < AHA T r NS # 10#2# UEIO# Z# 10. UKAREA# IER) 

C CALL USUCV < MIL AM "HR AD , H . WE I (J » NS , 1 # 4 ) 

CHAf«2)«*ll‘ 

CALL USUFH ( CHAR #12# AH AT r 10 . R'S # NS , 4 ) 

APPEND ZEROS TO FEEDBACK MATRIX AMD XFORM TO ORIOIMAL COL 
DO 701 X^l.NI 
DO 701 J»NOU# ORIOO 
F(I #.!) = FLOAT (0) 

READ<37fRL04>(<»<IIfIJ>,IJ=lfNT># IX-ltNS) 
READ<32'REC»S><<C<IX»IJ>f IJ»1 #N3> # II-l #N0> 

READ<32# REC=A) <<F(XI,IJ)#IJ“lfNS)fII**l#lll> 

CAL L USUFH < SHF - AFT # H # F # 10 # HI # NS # 4 ) 

Ctttttttt COMPUTE XFORM MATRIX M ft MINV ***************** 

DO 703 1=1 #NH 
DO 703 J s «l#NS 
MCI v.l)«FLOAT (0) 

IFC I .FO. J)M< X# J)**FLOAT ( 1 > 
ir(I«U:.NO)M(If J)=C(If J> 

703 CONTINUE 

CALL L 1 N V2F < M # NS # 1 0 , M I NV , IDO 1 . UKAREA r IF.R) 

C PAUSE 13 

C CALL UHUFH<SHM]NV|i*5tMINVflO#NStNS#4> 

C CALL USUrM<2HHf #2tl1#10fN3fNS#4> 

C ******** XFORM SYSTEM AND F MATRIX * ********************* 

CALL VMIILFE (hi A, NS » NS. NS» 10# 10# AHATf lOflf.R) 

CAL L VMIILFF <AHAT» MINIMIS i NS i IIS * 10*10»ALAM. 10# IER) 

CALL VHtILFFCMf D* NSf NStllli 10# 10# BLAH# 10 t IER) 

CALL VHIILFF - (CrHlNVtNOrNSt NSi 10# 10# CLAM# 10. IER) 

CALL VMULFFIF #HXtlV»MI #NS#NS# 10# 10# AHAT»10v IER) 

Dll H1H 1 = 1. NX 
DO 040 J=*l » NS 
FLAM ( I # . I ) « AHAT ( I » .1 ) 

040 CUHTIIHIE 

C********* DISPLAY XFORMED SYSTEM ********************** 

UI?nE(*»H01) 

001 FORMAT ( 1X« 3011 UISH TO DISPLAY TILDA SYSTEM?) 

kcadcs#*>kk 

I F ( KK . LE « 0 ) HO ft) 002 
WRITE! A# 704 ) 

704 FORMAT < IX • 2011 SYSTEM XFfJRIIED BY M#//> 

C PAUSE IS 

CALL IJSUFH (7IIA TILDA » 7 # AL AH # 1 0 # NS # NS # 4 ) 

CALL USUFH (7HD TILDA# 7# BLAH. 10# NS# HI #4) 

PAUSE 14 

CALL USUFH C7HC TILDA #7# CLAM# 10 » NO# NS i 4 > 

CALL USUFH <7HF TILDA. 7 1 AIIAT# 10»NI.NS#4) 

C ************** PARTITION THE SYSTEM ttttttttttttttttt 

002 DO 70S 1=1# HO 
DO 7015 .1=1 . NO 
ALAH1) (I# J)=ALAI1(I#.J) 

JF(I HTO.OT •HSMiiJ TO 023 
AL.AH21 ( I# J)»ALAH( I fHO#.l) 


701 CONTINUE 
C PAUSE 11 

C CALI. USUFH ( 1 01 IF APPENDED t r 1 0 , F # 1 0 . III . UR I HO , 4 ) 

CALL VHULFf ( F . VI) INV »NT % OR I HO . ORIOO .10.10. AM AT . 1 0 . X FR ) 
DCl 702 1 = 1. Ill 

do 702 j=i»orm*o 

F ( I # J ) “AIIA1 ( I # J ) 

702 CONTINUE 
PAUSE 12 

CALL USUFH (10HF XFORMEDI # 10#F# 10# Ml. OR ICO, 4) 

WRITE <32. RFC** A) ((K(IIrlJ)v X J«1 # ORIOO) #11 = 1. NX) 

C CALL USUFH ( 2IIF 1 . 2 # F # 10 # III # ORIOO # 4 ) 

C******** RETRIEVE URIHIIIAL SYSTF.M DATA ft STORE F ********* 
AW READ < 32 . REC=2 ) NS . MX » Ntl 

READ(32,REt =3)((A<II»I.I). 1.1=1 .NS) i XI»1 #MS) 

023 IF ( JIN0.01.NO >00 TO H24 
ALAH1 2 (I # .1 ) = Al.AM ( X » J I NO ) 

024 IF <( UNO .HT. NS). AND .(.UNO .01*. NS)) 00 TO 7015 


AL.AH22 < J » J ) “At AM <’X TNO. J4NU ) 

70S CONTINUE 
C PAUSE 40 

C CALL USUFH < /HAT 1 L 1 2 1 . 7 , AL AH 1 2 #10# NO . 1 2 # 4 > 

C CALL USUFH < 7IIATIL22 1,7# ALAM22 ,10,12,12,4) 

DO 70 A 1=1# NO 
HO 70A J=1,NI 
DLAM 1 ( I , J ) =DL AM < I » J ) 

ir < 1 1 NO . LE # NS > DLAM2 ( I , .1 ) »BLAH ( 1 1 NO , J > 

70 A CONTINUE 

DO 707 1=1. NO 
DO 707 JH.NII 
CLAMi ( I # J)=CLAH( I #.l) 

I F < J I NO . LE * NS ) CL AM2 ( I # .1 ) *CL AM < I , J I NO ) 

707 CONTINUE 

DO 7 OH 1=1 .NX 
DO 7011 .1*1# NO 
FTILKXi J)»AlfAT(IiJ) 

] F< J I NO. LE 4 NS )FT1L2< X #.J) “AIIAT ( I . .UNO) 

70S CONTINUE 

C******** TRANSPOSE AL.AH22 ft -At. AMI 2 TfJ ASS I UN OBSERVER DYNAMIC 
UR 1 1 E ( A # 022 ) NS # NO 

1122 FORMAT (IX, 3HNH« r 12 # 3X , 3HN0* .12) 

W0B5 =NS -NO 
DO 70V 1=1, NO 
DO 70V J» I# NODS 
ALAH12(I# J)»-ALAH12(X#J> 

70V CONTINUE 
C PAUSE 1A 

C CALL llSUFH<HK"ALAH12t #H,ALAH12. lOvNO.I'OBS. 4 ) 

C CALL IISUFM < 7IIALAH22 1 # 7# AI.AM22 ,10. NODS , NODS • 4 ) 

CALL TRANS! (ALAH22, NOBS# NOBS) 

C Al L TRANS 1 ( ALAM 1 2 # NO # NODS ) 

C PAUSE 17 

C CALL USUFH ( VII- ALAM 1 2Y I , V , AL AH 1 2 # 1 0 , NODS , NO , 4 ) 

C CALL USUFH ( HIIALAH22T I , H . ALAN22 #10. NODS # NOBS , 4 ) 

IF ( NO. LE. NOBS) HO TU 1131 
XFLF“1 



CALI. IJHWF'H ( VH- ALAN 1 2T I # V * ALAN 12 , 10 , NOBS .NO, 4) 
writf:<a,033>nobb 

HLU FORMAT ( IX# 21IJYOU MILS! SELECT WHICH. IS*. 42H OUTPUTS Wil l. BF IISFD Tl 
HEED THE OBSERVER! 

WR1TE(4,034> 

1134 FORMAT (1X#45HTHE OUTPUTS CORRESPOND TO COLUMNS IK -ALAM121) 

837 URITE<A.035> 

033 FORMAT ( IX, 40HSELECVED COLUMNS WILL FORM A MAI RIM THAT MUST BE HUH 
I INOtJI. AR ,//, 34HENTER OUTPUTS TO BE USED (INTEOEim 

RE AD ( 5 , * ) ( I LR ( I > . I * 1 . NOiiS ) 

no o:<4 j=«i .nobs 
no 03 6 i«i»nobb 

Al I AT ( I « J >«* AL AM 12 ( X . XLR(.H ) 

034 CONTINUE 

CAIL USWFM(4III< MODI .4. AHAT, 10. KOBO. KOBO. 4 > 

IER=0 

CALL LI NV2F C AHA T # NOBS # 1 0 # AL # I DOT . UK AREA * I ER ) 

If (1LR.EQ. 127)00 TO 037 
NOX*NOPS 
00 TO IMS 
031 NOX=NO 

no 030 I -1. NOBS 

no 030 J»1«N0 

Al IAT ( I . J ) « AL AM 1 2 ( I , J ) 

038 CUHIINUE 
045 COM TIKI IE 

C ******* STORE ALAM22 X AL AM 1 2 AO A fc » AND CHANCE NO ********* 

WR X TF: ( 20 . RL* l> I ) NOBS » NOX . NOBS , I POT i ZERO 

WRITE (20 1 RE02 ) ( ( AL AM72 ( XI . X J ) . I J =1 . NOBS ) . I X * i . NOBS ) 

URITE'20,RE(>3>((AHAT(II.XJ>, I.f»l #NOX> r IIM #KUBS> 

0*****< AiiSXON EIGENVALUES AND KIOENVEOlORS FOR OBSERVER ******** 

CL OSE (UNIT * HIT . STATUS" ' KEEP ' > 

MO 815 1=1. NO 
J=X F20 

CLOSE (UNIT- J) 

015 CONTINUE 

UR I T E ( 4 1 7 10) 

710 FORMAT ( IX. 20< lit*) t 32H ASSIGN El SEN VAUNTS FOR OBSERVER. 10 < lit* 
CALL HOPES 
WR I 1 L ( 4 . 7 1 1 ) 

711 FORMAT ( IX. 15 ( 111*) .3311 ASS t ON EIOEH VECTORS FOR OBSERVER# 14 (UN 
CALL M0HF3 

WRITE (4 .735) 

735 FORMAT ( IXi 33II0ISH TO REDUCE DA IN FOR OBSERVER?) 
l<EAIH5f*)KK 
IF (ICK «LE« 0)00 TO 734 
CALL MODE* 

734 CONTINUE 

C***t*t* RETRANSPOSE AI.AM22 X -ALAM12 X XPOSE F TO OF.T l. ********** 
NOX=ORIOO-NODS 
DU 712 1*1. NUBS 
DO 712 JolvffOX 
ALAM1 2(1 # J)**-ALAM1 2( I # J) 

712 CONTINUE 

CALL TRANS 1 ( AL AM22 . NS * NS ) 

CALL TRANS I ( AL AM 1 2 , NOBS . NOX ) 

C PAUSE 1 0 

C CALI. USWFM(7HALAM12I 1 7» ALAM12# 10. NOX. NOBS. 4 ) 

C CALL USUKM < 7HALAM22 1 # 7 # ALAM22 » 1 0 . NOBS , NUBS . 4 ) 


IUT-204NSII 

READ(IUT,Rt:i>4>(<F(II#IJ>#IJ=i#NS>#XI = l,NI> 

IF ( 1FLF.NE • 1 >00 TO 041 
no 03? 1*1 .NOX 
no 03? J* It NOUS 
1.(1# J) »EI.OAT Ml) 

03? CONTINUE 

PO 040 1=1, NS 
no 040 . 1 = 1 #NS 
L(lLR<I)..f)=F(I..ll 

040 CONTINUE 
CO TO 044 

041 HO 042 1 = 1 »NOX 
DO 042 J-i.NODS 
L( I # J)*F ( I r J) 

042 CONTINUE 
044 CONTINUE 

CLOSE ( UNIT » HIT . STATUS* ' KEEP ' ) 

CALL TRANS 1(1.# NOX . NS ) 

PAUSE 17 
IFLr*0 

CALL USUFI1 ( 2HL I # 2 » 1. 1 1 0 r NS i NOX * 4 ) 

C************* COMPUTE R *************************************** 

REAP ( 32 » REC«2 ) NS » HI . NO 
NOBS*MS-NO 

CALL VMULFF(FTXLi!.L#KI. NOBS# NO# 10. Id# AHAT# 10#XER) 

DO 713 1=1, NX 
MO 713 J»i .NO 

C R XS REPRESENTED DY KT1I.1 

F IX L 1 ( I , J > -FT1 L 1 ( I , .1 ) FAHAT ( I , J > 

713 CONTI Nt IE 


C ************** COMPUTE F **************** 


CALL VMULFf (LrALAMI 2. NOBS .NO. NOBS. 10, 10. AHAT . 10.XLR) 
| iO 714 1 = 1, NO PS 

no 714 j=i.Nnns 

FOBS ( X , J ) »ALAM22 < I # J ) -AHAT ( I r J) 

714 CONTINUE 

C*********** COMPUTE 0 *********************** 

CALL VMULFF (EOns, L, NOBS, NODS.NU. 10, 10. AIIAT , 10. XFR) 
CAt.l VMULFF <L , AL AMU # NUBS , NO , NO , 1 0 . 1 0 , AL AM22 ,10,1 ER ) 
1)0 715 1 = 1. NOBS 
PO 715 .1=1, NO 

C 0 XS REPRESENTED DY AI.AH22 

AL AM22 ( X , J ) *AL AM2 1 < I . J ) -AL AH22 ( X . J ) F AHAT ( I , .1 ) 

715 CONTINUE 

C *********** COMPUTE B- TILDA *********** 

CALL VHULF F ( L # BLAH 1 , NOBS , NO , HI , 10 , 1 0 . AHAT ,10, XER) 

DO 714 1=1, NOBS 
DO 714 J=1 » NX 

C AHAT REPRESENTS T-BTILDA 

AHAT ( I # .1 ) *■ BI.AH2 ( I # .1) -AHA T (X , .1 ) 

714 CONTINUE 

C********** display r.o.e.x tb ******************** 

CALL UBUF*M( VHMATRXX El #7#E0BS# 10, NODS, NOBS, 4 ) 



CALL USUFM ( 9NMATR I X Rl , ViFTILI » 10, MI ,NO, 4 > C 

PAUSE 20 

CALI. USUFM ( 9HMATR I X (M ,9, ALAH22,1<>,N0»S,N0,4> 

CAtL IJSUFIK 10HHATRIX T»l, 10, AHAT, 10, MOBS, til *4) 
r ******* CONSTRUCT TOTAL SYSTEM MATRIX ATILIJA ********* 

NS*NOBSTNO 

CALL VHIILFF(BLAM,FLAM»NS,NI,NS,10#10,AHAT,10, IER) 

|iO 710 I«1,NH 
no 710 J»1,NS 

AHA f ( I ,JI«AHAT< 1 1 J> 4AUiM< 1 1 J) 

71 H CONTINUE 

C CALL El ORF (AHAT , NS , 10 , 2 , UEIU , Z , 10, UiCAREA , IER ) 
f; CAL L USUC V ( 1 ONE 1 OENVALUE » 1 0 , UE 1 0 * NS » 1 # 4 ) 

CALL VNULFF<nLAM,Fm2,NS,NI,N0Bafi0#10#FTlU,10fIER) 
no 719 I“1 »NB 
LiO 719 .1-1, NOBS 
AHAT<1,.JTNS>"FYIL1<I,J> 

719 CONTINUE 

00 720 1*1 i NOUS 

no 720 j«i,nobs 

AHAT< X HIS ». IT NS) "HUBS ( I » J) 

720 .CONTINUE 

no 71 7 I -1, NOBS 
no 717 J«J,NS 
AHAT < ItllSf .D^KLOAT (0) 

717 CONTINUE 

I ROUGHS t NUBS 
CHAK’ ( 1 ) m ‘ AT1LUA TO* 

CHAR ( 2 > * • T AL ‘ 

C CALL USUFM i CHAR , 1 J , AHAT ,10*1 ROW , I ROU , A ) 

C CALL HI ORE ( AHAT , IKOUr 10,2,UEIO»Z, 10, UICAREA , I Ell ) 

C CALL USUCV< lOHEIOENVALOEr 10,UE10,NS, 1 ,4> 

C********* CONSTRUCT XFORH TO OET X IS U COORDINATES *** 

C FIRST MODIFY L TO BECOME M-|.li:i 
C 

MO 721 1*1, NOE'S 
IUI 721 J»lfNO 


XF < J*liT *NS )(»() TO 721 

XF < J, EO t X ) I. < X , J f NO ) =K LOAT ( 1 ) 

721 CONTINUE 
C 

C NOW CONSTRUCT M- INVERSE 
C 

NCI 722 1*1 , 1 ROU 
no 722 J-NST1 r IRON 
N1NV( X »J) "FLOAT (0) 

IF (X «K0« J) MTNVC X ,.l) "FLOAT <11 

722 CONTINUE 

HO 721 1*1, NOBS 
HO 724 J»1 ,NS 
MINV< I FNf»,.j) < c L ( I , J) 

724 CONTINUE 
C 

C CONSTRUCT M 


HO 725 1*1, IROU 
no 725 J-NSH , XKOU 
MCI * J)*FLOAT <0) 

X F < I . Efl . J > M < X , J ) *FI.OAT < 1 ) 

725 CONTINUE 

DO 724 I«1 iNOBS 
IUI 724 ,I<«1 , NS 

724. continue: 

CALL VMULFF (I. , M, NOBS, NS, NB, 10, 10,FT1L1 , 10, IER) 

00 727 I“i , NOBS 
HO 727 J-liNS 
NdTNftr.D-KTlLlCr, J) 

727 CONTINUE 
PAUSE 21 

C CALL USUFH C 2HM 1 . 2 , M , 1 0 , 1 ROU , I ROU , 4 > 

C CALL IJSUFHCtSHHTNVI ,S,H1NV» 1 0, IROW, IRON, 4) 

«mmi XFORH AHAT BY NEU M *************** 

CALL VMIII.KF (MINV, AHAT * I ROU r I ROU > IRUUv 10, 10,FTIL1 p 10, IER) 
CALL VMULFF (FT1L1 ,H, IROUt IRUU, IROU, 10, 10, AHAT, 10 , IER) 
C****t*4* CREATE NEU B-MATRIX FOR TOTAL SYSTEM ********* 
READCI2 ,REC*4 ) < ( H( 1 1 , 1.l) , I J=»l , N1 ) , 1 1 = 1 , NS ) 

BO 720 I- If NOBS 
liO 720 J*1,NS 

720 CONTINUE 

CALL VMULFF (L, ML AM, NOBS, NS • NI , 10, 10, F TIM , 10 » IER) 

NO 729 1*1, NOBS 
no 729 J*l, NJ 
BCITNRf J)"KTXLl(Xf J) 

729 CONTINUE 

C******* CONSTRUCT NEU C MATRIX FUR TOTAL SYSTEM ********** 
KEAD(32,REC*3) ( (C( XI r XJ) * XJ**1 , NS) , II"1 » NCI) 

DC) 7 .TO X"1,H0 

no /;ic» jtrHNSfiRou 

C ( I , J ) "FLOAT < 0 ) 

/.TO CONTINUE 
PAUSE 22 

CALL USUFH ( UHAHAT 1 » 5 1 AHAT ,10,1 RUU , IROU , 4 ) 

CALL EXORFCAHAT, IROUt 10,2, WEIS, Z. 10, WKAREA, IER) 

CALL (JSUCVC 10HEX SEN VALUE , lOiUETO, IROU, 1,4) 

CALL (TOUCH <1 ONE I SENVECTR, 10,7,10, IROU, IROU, 4) 

C****«*«* STORE NEW SYSTEM AHAT, fcr.C, NS, Nf , HO ************** 

UR I TE ( 20 , RLC= 1 ) IROU r N1 , NO , I DOT , ZERO 

WRIT! <20,REC~;i ) < <B(XI , IJ) • IJ=1 ,NI ) r II“1 , IROU) 

UR X TE < 20 , RE04 ) < (C(X 1, 1 J) , I J*1 , TRUU) , 1 1 - 1 , NO ) 

IUT=IR0UT20T1 

OPEN < UMX T » 1 LIT • F X LK« ' CURRNT ' , ACCESS 3 * DIRECT 9 ,RIXL =• 1 02 ) 

URX1EC XUf , REC~ti) < (AIIAT < II , X J) * 1,1-1 , IROU) , XX"i , XROU) 
C******** ENTER INITIAL CONDITIONS ANH XFORH TO FINN WCO) ******* 
URIT E<4 , 73 1 ) 

7T1 FORMAT < IX, 4511 ENTER INITIAL CONDITIONS FOR UR TO INAL STATES) 
REAB(5,# ) ( V0( I * 1 ) , I s *! ,NH) 
no S52 T»lTN(TBS 
no 052 . 1 * 1 , NS 
FTIL1 < X < J)w-MINV( ITNSf J) 


CD 



852 CONTINUE 

CALL VMULFF < F TILI . VO p NOUS p NS 1 1 p IO p J 0 . V p 10 p 1F.R ) 

no 732 X»1 pNUBS 

VO d 4 NB p 1 ) *V0 ( I f HO pI)-V(IpI) 

732 CONTI HOF 

,33 FUtm^faXtSlH use the following initial conditions in TIME RECF 
CALL 0HWFM(5HX«Dlr«»V0f lOtlllOUt It A) 

c ********* no TIMF RESPONSE *************************** 

WRITE(Ap4A3) 

4A3 FORMAT! IX'p25HW1SH TO DO TIME RESPONSE T ) 

READ(ti.*)KK 
IF <KK *I.K • OMNI TO 4A4 
CALI. MONK 4 
444 WRITE <4 » 734 ) 

734 FORMAT! 1Xp40HWISH TO REASSIGN OBSERVER AND TOY AGAIN? > 

RFAUC.if *)KK 
IF (KK.OT *0)00 10 AW 
Cttttt* RESTORE INITIAL SYSTEM INFO 
READ ! 32 • RFC-*) NS p N I p NO 

READ132tREC«3)<<A!ItJlf3*lfNSlf 1*1 'NS) 

READ ! 32 p RE<>4 ) < ! D < 1 p .1) » .1- 1 p NI ) p I* 1 » NS ) 

READ < 32 p REOti > < ! G < 1 1 J ) « J* 1 * NS ) 1 I * 1 1 NO ) 

WRITE! 20 r RFC* 1 DNSpNI »N(Jp IDGTp ZERO 
WRITE ! 20 r REC-2 ) < < A! I p J) p.l=l tNS) p 1*1 p NS) 
URITE!20pREt>3)((D(IpJ)pJ»lpHI>*I-ltNS) 

WRITE! 20pREC-4) ! !Cd p,I) pJ«1 pNS) p 1*1 pNO) 

RETURN 

end 

r;*t **************************************************************** 

SUBROUTINE TRANS! (At IMp IN) 

REAL A! 10 p 10) r AT < 10» 10) 

DO 10 X=1pIN 

no io j* 1 pin 

AT! Jp I )»A( I p J) 

10 CONTINUE 

DO 20 I«1 r IN 
NO 20 J=IpIM 
A< X ? J> C AT< I p J) 

20 CONTINUE 
' ^RETORN 
END 

c****************************************************************** 

SUBROUTINE VAST 
INTEGER OR I 00 

REAL VO(iOpiO)pl-RORO!iO)pLIOrNi(10).VFS<10p 10) » VDES! 10. 10) 

1 r ALPHA I ! IO p 10) pREL! 10p 10 ) p ALFHAR !10p 10) tKIODir! JO pIO) t 
2VJ *< 10p 10) t BLAN2! 10t 10) iR< 10 p 10) » VOR! 10* 1 ) r IHL< 10 p 10) 

3pMI ! 10p 10) p VOI ( 10p 1 ) pVPSTMP! 10.1) ? VEST PC 10 p 1 ) p AL.AM22! IOp 10) 
RIAL W ! 1 0 1 10) pM ! IOp 10) i VINO! 10 p 1(0 pF < 10# 10 ) p Alt AT ! 10» 10) 

REAl. XX ! 10. 10 ) p VA < 20 ) p E C 20 ) p X < *0 ) . I.RF ( 10 ) . L lii< 10 ) r UJUO ) 

REAL DI0!20p20)«ni0IN0!20>20)pWK11I0<4.;0) 

REAL A<10» !<>).»< 10« 10) pC! 10. 1(0 . WLAREA! 130 ) 

COMMGN/B Y8/A p D p C p ZERO . I DGT » NS p Ml p HO 
CUMMUN/AUO/F . Al I AT /E I S/LRE p 1.1 M/PAR/ At /GR/G 
('OMMGN/VEC/VA*E* X*UJf U*XX p V tVINV 


CUMMOM/RU/UR T SO . VO p VFS p VDEB p BI.AM2 p LRClRU p L 1 URS p AI.AM2 
I2=()RIG0-NS 
10=0 

i 111=1041 

IF ! ADSO.IM! 10) ) .ST •ZERO) ICMFLX=i 
C **************** CREATE LAMBDA'S I ******************* 

DO 10 1=1.12 
no to j«i*X 2 
REL! Ip J) -FLOAT (0) 

XFO.FM. J) REI.d* J)=LREdO) 

10 CONTINUE 

C CALL USWFM ( 3HREI. p 3 p REL p 10 p 1 2 p 12 • 4 ) 

C**************** CREATE LAMBDA- 2 *********************** 

HO 20 1=1 p 12 
HO 20 J= 1 r 1 2 
ALFHAR ( 1 1 J ) »AI. AH22 < I p .1 ) 

20 CONTINUE 

C CAt.L USWFM ( AHALFHAR p A p ALPHAS p10pI2»I2p4) 

C **************** REL -Al.FHAR ********************* 

DO 40 1=1.12 
no 40 J*lfI2 

EIGDIFUt J)*RELdf.))-ALPHARd».t) 

40 CONTINUE 

C CALI. UBUFM! AHEIGDIFp ApEIGDIFp IO p I2 p 12.4 ) 

C PAUSE '40' 

IF!) CMFLX • NE • 1 ) 00 TO 34 

C************* CREATE UKtl ******************** 


00 71 I"1 .12 
DO 71 J*1#I2 
IML < I . J ) «FI.OAT ( 0 ) 

IF! J .E0« J) IML! I pJ)*LJM< 10) 

71 CONTINUE 

C************* CREATE BIS ********************** 

BO 72 1*1 » 12 
L>0 72 .1*1 1 12 
Diodp.n-EisniFdpJ) 

BIGd4I2pJ)=-IHLd.J> 

BIG! I p J4I2) »IML d p J) 

B1GU4I2. J4X2>=El0BXFdp.» 

72 CONTINUE 
122 = 2*12 

CALL LINV2F < BIGr 122 p20p BIGINV* I DO l pWKMGp 1ER) 

DO 73 1=1 p 12 
DO 73 J=i . 12 
AL.FHAR (I r J ) »BI G I MO ( I * .1 ) 

IML< I p J)=HIGINVd p.1412) 

73 CONTINUE 
GO TO 35 

C*********** TAKE INVERSE OF EISDIF ************* 

34 CONTINUE 

CALL LINV2F!EISDIF» 12* IOp ALFHAR p XDST . WK AREA t X Eft) 
f;** ******** FREMULT BY V12 ************************ 

35 CONTINUE 

BO 41 1=1 »NB 
HO 41 J-1.X2 
VI 2! 1 p.J)>*V()d p J4NS) 

41 CONTINUE 
C PAUSE ' 4 1 ' 



C CALL USUE ri < 71 II N VERSE 1 7 . ALPHAR 1 10 . i 2* 12 , 4 ) 

C Al L VHUI.FI ' ( V 1 2 r ALPHAR * NS . 1 2 * 1 2 , 1 <» . A 0 * REL . 1 0 . XEK* ) 

C PAUSE 'REL' 

C4 ********* P0S1MUI.T BY BLAH 2 ******************* 

CALL VHULFF (REL. BLAH?* HU. 12. N1 • 10* 10, ALPHAR? 10. IER) 
f; PAUSE 'MAH2' 

i : ************ POSTMHI T BY F ********************* 

CALL VMIILFF < Al PHAR . F t NS * M X t NS * 1 0 » 10 » PEL » i 0 . X Eli ) 

0 PAUSE 'F' 

C************tAH RESULT TO VI 1 
BO SO 1=1* NS 
BO HO J“1 *NS 
R < I r J ) =KEL ( I . .1 ) 4 VO < I * J ) 


HO CON I I HUE 

C********** POST MULTIPLY BY ASSIHNEH EIGENVECTOR ********** 

BO AO 1=1. NS 
VQRU. 1>«»V(I*IU> 

AO CONTINUE 

C CALI IISUKM(3HVQR.:SiV(IR.10»NS*t.4) 

C CAIL IISUF‘M( 1HR* l.R* 10*NS.NS*4) 

CALL VHULFF C R * VOR • N8 » NS • I f 10 1 Id • VFB1 HP « 1 0 1 XER I 

BO 70 1*1 .NS 

VKS( I » XU)=VE8THP( I . 1 ) 

70 CClNIlNUE 

XF(ICMPI.X.flE.i) 00 TO 111 

C*********** PREMOLT BY VI2 ************************ 

CALL VHULFF ( VI 2 » IMl. . NS . 1 2 . X 2 . 1 0 * 1 0 » REL » 1 0 » I ER ) 

C*«« *********** POBTMUl.f BY BLAM2 ****************** 

CALL VMOI Ff < REL . Bl. AH2 . NS * 1 2 . N1 ► 1 0 1 1 0 . AL PHAI .10. IER) 
C**************** POSTMULT BY F ******************** 

CALL VHULFF ( AL PHA X » F» NS f N1 * NS . 1 0 . 1 0 * MX 1 1 0 . XER ) 
r; ************** COMPUTE FZIOENVECTDRK ************** 
no HO 1*1 t NS 
VQI ( X* 1 ) n V( I . XU41 ) 

HO CONTINUE 

CALL VHULFF (Hit VOX . NS.NSt 1* 10 1 lOvVFHTHPv 10 » 1F.R) 
no 90 1*1 .NS 

VF6 < X » 10) «VES ( 1 1 XU ) -VEST HP ( I » 1 ) 

90 CONTINUE 

CALL VHULFF ( R . VO I » NS . NS » 1 . 10 , 1 0 * VF6THP 1 10» IER) 

CALL VHULFF (MI » VOR.NH.NB. I * 10* 10. VFS1P* 10. IER) 

BO 110 1*1 * NS 

VES < I * IU+ 1 ) « VFBTHP < X . 1)4 VK STP (1*1) 

110 CONTINUE 

ia=um 

111 CONTINUE 
ICHPLX-O 

IF < KM. T. NS) SO TO 1 

RE TURN 

ENB 

C ***************************************************************i 

C*********************** ****************************************> 

SUBROUTINE ROCOST(CJ) 

INTEGER OR I SO 

REAL XX(10*10)*VA(20)*F:(20)»X(20)*LRE<10)*I.IM<10) iUJ< 10) 
REAL U( 10* 10). VC 10 *10) .VXNVU0.10) t AL( 10. 10) . AI.AM22( 10. 10) 
REAL VES ( 10* 10) *VI)ES( 10* 10) • BI.AM2( 10* 10) * VO ( 10* 10) 


REAL A(lO*iO).B(lO,iO)*caO,lOr;rKim7ja0).LUJRO(iO) 

CONMON/S YS/A . B * C . ZERO * I DO f » NS • N X » NO 
COMHON/VEC/VA . FI * X * U J * U * XX . V . VX NV 

LOMMON/RO/OR X HO * VO * VFS * VBES * BI..AH2 * LRURO * L X ORS . AL AM22 
COMMON/I: X 0/LK‘E * L X M/IJAR/AI. 

ctt ******* calculate actual V- xn FOI L SYSTEM ******************** 
CALL VACT 

C CALL IIBUEM(3HVF'H. 3. VF’S* 10*NS*NS» 4 ) 

C CALL USUF* H ( 4I1VDFZS * 4 * VBES » 10 * NS * NS > 4 ) 

c********* calculatf: cost function ****************************** 

CJ=0 

J=0 

5 J*J41 

1*0 

10 1*141 

C JTF.MP“ ( VES < X * J ) - VBES ( I * J ) ) * *2 
(;JCJTFMP*AL(X» J)4CJ 
XF(ABS(LXM(J)).LE.7ER0) HO TO 20 
C J=C J4L JTt MP*AL ( I * J ) 

CJ=CJ4< Vl*8< I » J41 )-VDFZS( I * J41 ) )**2*2*AI.( X * JF1 ) 

XF < X * NE «NS) HO TO 10 

.Kill 

20 if(i.nf:*ns) ho ro io 

IF(J.NE.NS) HO TO 3 

RETURN 

ENB 

C**************************^*************** ********* ***************: 
C********************************* ******************* **************; 

SUBROUTINE ItOHRAD 
INTLliER ORXHO 

RF AL AUX1 ( 10 . 10 ) . AIJX2 (10*10). AIDS! (10*10)* PVCRE (10.1). PVO X ( 1 0 * 1 
1VAUX1 ( 10* 1 ) * VAUX2( 10* 1 ) »BETAR( 10*1)* L'ETARi 10*1)* 

2V11 ( 10.10)* BETA I (10*1) .ZETAI (10* 1 > . RIU 10, 1 ) *WK 10* 1 ) 

J » RO ( 1 0 * 1 ) * OR ( 1 0 * 1 ) * 0 ( 1 0 * 1 0 ) » Al. ( 1 0 * 1 0 ) 

REAL 00(10*10), I.RORO (10), I.XURO (10), VFS (10*10), VBES (10,10) 

1 , ALPHA! ( 10* 10) , Rt l.( 10, 10) , ALPHAIU 10* 10) tKXUDXF ( 10 , 10) . 
2V12(10*10)*BLAH2(10*10) ,R( 10* 10) »VOR( 10* 1 ) • 1ML ( 10* 10) 

3,MX( 10*10) * VOX (10,1), VESTMP( 10*1), VKSTP( 10, 1) 

REAL ML ( 1 0 * 1 0 ) , III. ( 10 * 1 (» • NL C ( 1 0 * 20 ) * Pl.C (10,20), MI C (10,20) 

REAL STAR (20* 20 ) * NL ( 1 0 , 20 ) * RL. ( 10 , 20 ) 

REAL U<10,10),V(10»10) *VINV(10> 10) *F( 10* 10) * AII.VT ( 10 * 10 ) 

REAL XX ( 10* 10) » VA(20) *E(20) * X(20) *I.RE( 10) *LXM< 10) ,UJ( 10) 

REAL A(10*10),B(10, 10 ) i C ( 10 , 1 0 ) » UKAK'EA ( 1 30 ) 

REAL If X U ( 20 * 20 ) * B 10 1 HV ( 20 » 20 > * UKIUS ( 4*0 > * At. AM 2 2 (10,10) 

COMMON/ S YS/A »B,C» ZERO * I HOT * NS , NX , NO 
COMMOII/NSPA/ML , NL , NLC * Pl.C i ML C * STAR , RL , RL 
CUMM0N/AUX/AUX1 , A0X2* AlJXit 
COMMON/ At JH/F , AHA T/E X H/LRE , 1. 1 M/PAR/AL /tiR/ft 
COMHON/VEC/VA * E« X * WJ* U> XX. V* VI NV 

CUMMON/RO/OR XOII * VO * VES * VBES * BL All 2 . I.RORO * LI OKU , ALAH22 

JZ=ORXSO-NS 

1FLAH-0 

ns ioo x«i *ns 
no too j«ifKi 
C(l» J)«FI.0Af(0) 

100 CONTINUE 
1 i«=o 

110 11=1141 



10*0 

120 XU=IUfi 
10*10120 

iF<Amui.nuiun*ar,zt'Rcnm» to ;m 

READ < I U » REC<I > < ( NL < I*. .1 ) » .1* 1 1 Ml ) * 1 « 1 1 NS ) 
REAn<XU.REC-4)(<HL(I..I).J-l.NI).I=l.NI) 
oo rci 35 

34 IS-NS+NX 
NI2-2*N1 
N32=NS*2 
INS-NSfl 

read< iu.rei>3) < <nlc< i. . j»i , is> » i«i .ns) 

READ<IU.RE(>4>< <PLC< I » .1) . ,1=1 . IS) • Z-l . NS* 
REAVuxu.REc-5><(hu:a..)>»j-i.is>,i=i.Kx> 

Rt An<rU.REC‘'A><<m.(I..J).J-l.NI2>.I‘»l,Kft> 
ri.au ( i u » mri:*7 > < < rl c i » j> » J- i . N 12 ) » i« 1 . no > 

35 continue 

C PAUSE 'PFX' 

FALL PFXUX.IU.XFI.AU) 

C PAUSE * PFX * 

CALL PVP(lI.Itl) 

C PAUSE 'PVP' 

If <AB5<LlM<IU>>.LE.ZERU)til) TO 12s 
XFI.AU2«1 
IFLACI =1 
125 CONTINUE 

C************** CAICUI.ATK ALPHAR *********************************; 
{ ;**************«* CREATE LAMBDA -II I ******************* 

Pit 10 X-1.I2 

no 10 ,1-1.12 

REl.C I » .1) “FLOAT CO) 


IF < I * EO * J ) REL < I * J >=LRE ( 10 > 
tO CONTINUE 

C**************** CREATE LAMBDA'* £ ********************* <* 

DU 20 I«i .12 
DO 20 .1—1 r 1 2 
ALPHAR ( X * J ) =Al.AH2? (If J) 

20 CONTINUE 

C**************** REL*' ALPHAR ********************* 

DO 40 1=1 r 12 

no 40 j-i.xm 

EIGDIFC I • .!)=RLL< I . .f)~ALPHAK( I » .1) 

40 CONTINUE 

irrarLAOi.Ni:.i) co to t;i4 
c ************ create Lint I mmmmmi* 
no n 1 - 1.12 
no 7i j-i.ia 
imlcx. j)«n.0Arco> 
ir( /arc. j)im.u» j)=lim( iu) 

71 CONTINUE 

c *************** CREATE BID ***************** 

1)0 72 1 = 1*12 
DO 72 .J-1.X2 
mu< i * j)=EiuniF< i » j> 

DIG<H 12. J> — IMLCX. J) 

IUG<I»JtI2)-IKL Cl, J) 

0IGmi2. JTI2)'«El(iDIF<I,J) 


72 CONTINUE 
122-2*12 

CALL LINV2F(niUr 122 . 20 . B10INV. I NOT . UKDIU. IKR) 

DO 73 1=1*12 
DO 73 .1-1*12 
ALPIIARI I* .J)«BIOINV<Ir J) 

imlu.j>«bio<i.ju2) 

711 CONTINUE 
CO TO 136 

C*********** TAKE INVERSE UK EICDIK ************* 

134 CONTINUE 

CALL 1. 1 N V2F ( E I CD I E * 1 2 * 1 0 » ALPHAR t 1 HOI » UKAREA * I KR > 

C ********** PREMULT BY V12 ************************ 

136 CONTINUE 

BO 41 1-1 .NS 
BO 41 .1-1.12 
V12CI..I)*=V0(If JINS) 

41 CONTINUE 

CALL VMULFFCV12.ALPIIAR.NSt 12* 12. 10 1 10. REI. .10. XER) 

C ********** POSTMUl.Y BY BLAH? ******************* 

C A l L VMIII FF < REL » BL AM2 . NS • 1 2 1 Ml • 1 0 . 1 0 1 ALPHAR . 1 0 . X ER ) 
f; ************* CALCULATE NETAR ********************************* 
DO 130 X-l.NS 
PVORE (I r 1 > wAUXl C I r II) ) 

130 CONTINUE 

CALL VKIJLFF < F . PVPRE # HI * NS . 1 « 1 0 » 1 0 . VAU XI. 10. 1 ER > 

DO 135 1=1. NS 
VUR < I » 1 ) ** V (1.1(1) 

133 CONTINUE 

CALI. VII Ul. FF ( U . VRR . N X » NS » 1 . 1 0 1 1 0 . VAUX2 . 1 0 . I ER ) 

HO 140 1=1. NS 

BETAR ( X . 1 > »‘VAUX 1(1.1) TVAIJX2 < I » 1 ) 

140 CONTINUE 

C************* CALCULATE 2ETAR ********************************4 
BO 145 1=1 .NS 
1)0 143 J=l»Nfi 
V1KI. J>-VO(I.J> 

143 CONTINUE 


4 ZE1AR, 10. XER) 


no 150 1^1 .ns 
RR (1 . 1) =ZE TAR ( 1 . 1 ) 1 VAUX 1(1.1) 
130 CONTINUE 


-- IMI Ijr AO/ 

C**t ************$ CALCULATE ALPHAI *********************&tttlt^± 

«'«««»' PKEH.u r Br via 

..... ,:AU - W1MUTfVia.IHL,M8,12.ia,l0.10.REL.I0/XER) 

CUKIHMttttl POSlMlIkT BY BI.AIi;! HtltlttllHItlHt 

..... ,:AU - VMMUT<REL»BI.AMa.M«.Xa,KX,10.10.ALPnAI,lV,n-(0 

Ct*t**t**$***t**t* CAkCDkATK BET/iI ttttttttitttttintt ***»*«*« 

DO 155. 1=1 .NS 


155 


PVIII ( I . 1 ) ~AUX1 ( J . 1011 ) 

CONTINUE 

CALL VMIII. FF < F . PVR I . NI . NS . 1 . 10 » 1 0 » VAUX 1 . 1 0 . ICR > 
DO 160 I = 1 > NS 
VUIU. 1)«0(I, mil ) 



160 CONTINUE 

CALL VMULFF <U* VOI *NI »NS * 1 * 10* 10* VAUX2 * 10* IER) 

1)0 16!5 1*1 *NS 

BETA X < I r 1 ) * VAUX 1 (1*1) +VAIJX2 (1*1) 

UH CONTINUE 

c***tt*t***ttt compute zetai 

CALL VMI JUT ( V 1 1 • PVO I * NS * NS • 1 * 1 0 * 1 0 , ZETA I » 1 0 * X ER ) 
C****M**#t*lt** CALCULATE 00 ***tt**ttt*«t*t*t***tt*t***t* 

CALL VMULFFC ALPHA I * BETAI *NSi NX * l*10i 10*00*10* IER) 
r;***ttt*t****t*t CALCULATE OR *t*ttt*t*tt**t**t**tt**t*tt** 

CALL VHIfLFFCALPHARr DETAI *NS*NI * 1 * 10« 10* HR* 10# XER) 

1)0 170 X»1*NS 

0R( I * 1 ) "ZETAI <1*1 )4RR( I * 1 ) 

170 CONTINUE 

C ********** SEE IF THIS IS SECOND TIME THROMOH FUR COMPLEX tilt 
XI (IFLA02*NE.1> 00 TO 1 AH 
1A7 JJ»XO 

00 TO 16V 
l&ll JJ* 1011 
1AV CONTINUE 

C***t***ttt* CONFUTE GRADIENT FOR -REAL CASE mtittttttttttttttt 

ir-o 

17S IP*IP11 

G( X I * JJ) *0< IX * JJ)+(UFH( IP* 10) -ODES ( IP* 10) X <RR< IP* 1 ) 

1*2*AL (IP* X(l) ) 

IF(IP.NE.NS) CO TO 175 
XF( IFLAOl «NE« 1 >01) TO 203 

C**t**ttt*« COMPUTE GRADIENT FOR COMPLEX CASE ttmttttttttttttl 
1P*0 

100 IP*IPU 

om*jj>-o(iif.;.;>-(VFH(xPfXO)-v/)E«(XP*xf/))»<fiRaPri) 

1 1 2 * AL (IP* 10) )4( ( VFIi ( IP * 10 * 1 ) -ODES ( IP * 1(14 1 ) ) * ( OR < I P * 1 ) 

2fK0( IP* 1 ) )*2*AL( IP* I0F1 ) ) 

IF(XP.NE.NS) 00 TO 100 

C********** CHECK IF THIS IS FIRS! TIME THROUGH ***************1 
1F( in.A02.NE • 1 ) 00 10 201 
JTLA02~0 
DO 105 1*1 * NS 
DO 105 J ! « 1 * NS 
AUX1(X*.I>«AUX2(X*J) 
lUti CONTINUE 

DO 1V0 I«1*NX 
DO 1 VO 


U ( 1 » J ) - Al 1X3 ( X * ,| ) 

1V0 CGNUNUE 
00 TO 125 
20 J CONTINUE 
10=1011 
203 CONTINUE 

ir(iu.i.r.NS)00 to 120 
II (U.LT.NI)OO TO 110 
CALL PBNURM(NXiNS) 
RETURN 
LND 



124 


APPENDIX B: DESIGN SESSION EXAMPLE 


********************************************************************** 
********************* SPECTRA!. ASSTKNMENT PACKAttK ******************** 

ENTER DESIRED NODE OF OPERATION > M()HE«0 » i > 2 > . . . r 9 5 
? 1 

************************** NODE i {DATA ENTRY **********,*************** 
**********ENTER OR CHANCE SYSTEM PARAMETERS! 


PREVIOUS VALUES? 
? 1 


n«= 7 

NX* 2 

NO* 4 3 DOT- 9 

ZERO- . 

oooooo i ooooo 

W I till 
0 

MATH 

TO CHANGE? 





XX A l 

X 

6 

7 

:i 

4 

tv 


••.1540001- +00 
--•320000C-01 

♦154000E+01 

•OOOOOOE+OO 

*• . 420000+1-02 

♦OOOOOOE+OO 

- ♦ 744000E+00 

2 

-.VVAOOOE+OO 
« OOOOOOE+OO 

-.U7000E+00 

♦OOOOOOE+OO 

-.295000E-03 

•38AOOOE-OX 

♦200000E-01 


♦ 24V000F+00 

- ♦ j. j 2000 c i oa 

-♦ 520000E+01 
.OOOOOOE+OO 

-. iOOOOOE+Oi 

♦OOOOOOE+OO 

• X? 7000+“ +00 

4 

•OOOOOOE+OO 
♦ OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

♦ ioooooE+oi 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

tv 

♦OOOOOOE+OO 
• OOOOOOE+OO 

•*00000 OE +00 
♦ OOOOOOE+OO 

.OOOOOOE+OO 

♦OOOOOOE+OO 

- ♦ 200000E+ 02 

6 

♦000000E100 
-.25000011 102 

.OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

. OOOOOOE+OO 

♦ 0000 OOE + 00 

7 

ursH 

• 500000EJ 00 
« OOOOOOE 1 00 
TO CHANGE? 

.OOOOOOE+OO 
- . 500000E+00 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 
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MATRIX » i 



I 

2 




A 

* OOOOO OE f 00 

♦OOOOOOE I 00 




•) 

« OOOOOOE f 00 

♦OOOOOOETOO 




S 

•OOOOOOETOO 

•OOOOOOETOO 




4 

• OOOOOOE tOO 

•OOOOOOE tOO 




tv 

. 200000 ET 02 

« OOOOOOETOO 




A 

« OOOOOOL f 00 

•250000E+02 




7 

WISH 

► ('i 

4 OOOOOOE ! 00 
TO CHANGE? 

♦OOOOOOETOO 




MATRIX C J 

.1 

A 

2 

7 

;j 

•1 

!5 

1 

♦ AOOOOOETOA 
•OOOOOOE fOO 

♦OOOOOOETOO 
- ♦ 100000E+01 

•OOOOOOETOO 

♦OOOOOOE tOO 

♦OOOOOOETOO 

2 

* OOOOOOEt 00 
« OOOOOOETOO 

•OOOOOOETOO 

•OOOOOOETOO 

♦ lOOOOOETOi 

•OOOOOOE tOO 

♦OOOOOOETOO 

;? 

• OOOOOOE TOO 
♦OOOOOOETOO 

• AOOOOOETOA 
♦OOOOOOETOO 

•OOOOOOE tOO 

♦OOOOOOETOO 

♦OOOOOOETOO 

4 

WISH 

♦ OOOOOOETOO 
•OOOOOOETOO 
TO CHANGE? 

•OOOOOOETOO 

♦OOOOOOETOO 

•OOOOOOETOO 

♦lOOOOOEtOl 

♦OOOOOOETOO 


? 0 

WISH TO EXIT FROM THIS MODE? 
r A 

4c**^******^4c***^c*^*)jc«*«***< EXITING MODE I ******«***#*#*'4f>J:*i)c***^4:^* 
TERMINATE THIS RUN OR SELECT NEXT MOPE* 

WISH TO TERMINATE? 


^***:U;U:m:m***** SPECTRAL ASSIGNMENT PACKAGE *:m:m****aom*;r*** 

ENTER DESIRED MODE OF OPERATION, MfJ))E=0, i » 2 , * . * r 9 1 
? 9 

WANT TO ENTER NEW ORIGINAL EIGENVECTORS? 

? 0 

v matrix: 



1 

2 

3 

4 

t7 


A 

7 




A 

• SGS02AET00 

-• 1 19002ET00 

-♦8tvl902E~0i 

♦2744S1ET00 

♦37417AE-01 


♦ J 2V240E-02 

♦OOOOOOETOO 





♦ 1 1 i 44 SET 00 

♦ 309A71E TOO 

.A25909E--0S 

.25A35 8E-0I 

-«3AA4A5E-05 


♦ 5tv:iC02E-04 

•OOOOOOETOO 




s 

-•♦A03S27ET0A 

-.ssyysan+oo 

♦405A93Ei01 

— • AtvtvA 1 2E "0 1 

.t svvou;~o:i 


•4AAA52E-01 

•OOOOOOE too 




4 

-,22syoyE+oo 

♦G2VA27ETOO 

-.373754E TO! 

.7iS357Et01 

. 89VH02i:-0:i 


-•• JBAAAiE-02 

•OOOOOOETOO 
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U , OOOOOOt. + OO 

♦OOOOOOE+OO 

♦ OOOOOOE+OO 
• OOOOOOE+OO 

« OOOOOOE + OO 

♦ooooooe+oo 

♦100000E+0X 

, ooooooe+oo 
♦ioooooe+ox 

.000000+1+00 
, OOOOOOE+OO 

1 

♦OOOOOOE+OO 

.OOOOOOE+OO 

♦OOOOOOE+OO 

7 , 18AV97E-02 

~,26375SE“04 
WISH TO CHANGE ANY 
0 

alam; 

1 

/> 

“* 1 150215 ill +00 
♦ joooooE+o.t 
VALUES or V? 

• 727SS6l>01 

. 2 796 06 E +00 

- . 959426E-03 

y 

7 

a 

4 

5 

X -.001905 E-01 

• OOOOOOtriOO 

♦12A948E+01 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

2 • . i26948E+01 

.00000011100 

-•881901E-01 
, OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

;? .ooooooe+oo 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

~,X0SS45t: + 0i 

.OOOOOOE+OO 

♦OOOOOOE+OO 

4 ,000 oooe+oo 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

♦OOOOOOE+OO 

-.91A482E-02 

♦OOOOOOE+OO 

sS « OOOOOOE+OO 

, ooooooe+oo 

.OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

- .200000E+02 

0 .OOOOOOE+OO 

-•♦ 2 S 0000 E +02 

.OOOOOOE+OO 

♦OOOOOOE+OO 

♦OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

7 .ooooooe+oo 

♦ooooooe+oo 

WISH 10 DISPLAY RCl 

1 :i 

♦OOOOOOE+OO 
- , SOOOOOE+OO 
MODEL.? 

♦OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

REDUCED ORDER MODEL 




MATRIX a; 

X 

n 

s 

4 


x -.amvoiE-oi 

•126948E+01 

.OOOOOOE+OO 

♦OOOOOOE+OO 


2 » • J2A948E+01 

- • SSS 901E-01 

.OOOOOOE+OO 

.OOOOOOE+OO 


;5 .oooooon+oo 

.OOOOOOE+OO 

-♦10SS4SE+01 

. OOOOOOE+OO 


4 . ooooooe+oo 

MATRIX Di 

1 

.OOOOOOE+OO 

♦OOOOOOE+OO 

-.91 64 82E -02 


x ■■.iAai > ;5SE+oi 

1 A4440E-01 




2 ,A5oy2;5t:+oo 

,7229 09E-01 




J -*.27/2i5E+00 

- • 28 780 7E + 0 0 




4 ~.2A975;5E-I00 

MATRIX Oi 

1 

J5275IF+00 

2 

:s 

4 


1 .385 115J.E + 00 

.31249iE*0i 

-.I5794AE+00 

-.512K08E-02 
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:> - * xo3327E+oi «.;jsy 952t:+oo ♦ 1056 ? 3t:+ox - . asmaxze - ox 

3 . iii64uu-oo *3<>?67 xe+oo *62svo?e-ox *256358E--oi 

4 • . 225909E+00 * 829627E+00 - * 373754E+0X ♦ 7X5357E C 0i_ 

:tt:m****:m****$* REDUCED ORDER EIGENVALUE A3SXGNMi'NT***:ME*4*S* 
^^*^^^*^^^**1:*#*^^** mde ^eigenvalue assignment zx%***%*%*%*%% 

********** EKTIiR nit CHANGE EIGENVALUES! 


PREVIOUS VALUES? 
r 0 

LAMBDA Xi 

REAL= * OOOOOOE+OO 1 MAG* * OOOOOOE+OO 

WISH TO CHANGE? 
r .1 

ENTER NEW VALUE(S) * 

r -i.s i * 5 

LAMBDA 2 i HEAL* - . X50000E+01 t I MAG* - * 150000E+0X 

NEXT eigenvalue: 

PREVIOUS VALUES? 

? 0 

lambda ;*: 

REAL* * OOOOOOG + OO I MAG* *000000E + 0<> 

WISH TO CHANGE? 

? X 

ENTER NEW VALUE (S ) i 
? -*2 i 

LAMBDA 4 1 REAL* -.200000E+0X rXMAG* -• XOOOOOE+OX 

WISH TO EXIT FROM THIS MODE? 

? i 

*************************** EXITING MODE 2 ******************* 

ENTER DESIRED PARTIAL EIGENVECTOR ASSIGNMENT 
EIGENVECTOR V X 
r 20 0 A 7 0 0 0 0 
EIGENVECTOR V 3 
? 0 0 0 0 20 0 -0 -4 

USE THE FOLLOWING V MATRIX FOR INITIAL ASSIGNMENT 
REMEMBER WHICH V ARE COMPLEX 
INITIAL. GUESS FOR Vi 

X 2 3 4 

X . 406903E+02 . 695530E+0X - . 305037E LOO . 343961 E +00 

2 - ♦ X36521 E+OX , X94037E+02 ;■ • 90005XE+00 - . 969533E- OX 

3 , S24X2XE+02 * 34901 5E+0X *47089511+01 * 70922XE-0X 

4 *0i 0096E+01 ~ « 207752E+00 * X48707Ef0X - * 500000E + 00 

***********4 c******** MODE 3 i EIGENVECTOR ASSIGNMENT ************: 

********** ENTER OR CHANGE EIGENVECTORS i 


PREVIOUS VALUES? 

? 0 

EIGENVECTOR V X i (REAL.) ( I MAG ) 

, 0 0 0 0 0 0 E + 0 0 * 0 0 0 0 0 0 1 : + 0 0 

. OOOOOOE+OO ♦ OOOOOOE +00 

* 0 0 0 0 0 0 E + 00 , 0 0 0 0 0 OE f 0 0 

. OOOOOOE+OO ♦ OOOOOOE+OO 

WISH TO CHANGE? 
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7 l 

ENTER A NEW DESIRED VECTOR i 
? 40.7 6,96 -X.37 19*4 .12,4 3*49 0.10 -*.?00 
COMPLEX VDJ 


* 4 8 7000 E +02 
. 194000E+02 
ACTUAL VECTOR i 


--♦i;?7oooE+oi 

.349000E+01 


• .t 24000E+02 
- • 200000E tOO 


.BiOOOOE+Ol 


♦402275E+02 
•1V8133E+0? 
ERROR VECTOR i 


♦ 12) 094E+01 

•373668E+01 


♦127B16E+02 

-•&37207E-02 


•ooci30E+oi 


♦472546E+00 

-*4X3257EtC0 


-•2S8OV4E+01 

-•246684E+00 


LENOTH OF THE DESIRED VECTOR « 
LENGTH OF THE PROJECTED VECTOR 8 
OF THE ERROR VECTOR 
ERROR ACCEPTABLE? 


LENOTH 
10 THE 
? J 

EIGENVECTOR V 21 


NEXT EIGENVECTOR) 
EIGENVECTOR V 3) 


(REAL) 
•4B7000E+02 
-♦137000E+01 
• 1 24000E+02 

* BiOOOOE + Oi 

(REAL) 

* 00 000 OE +00 

* 00 00 00 E +00 
♦ 00000 0E+ 00 
♦OOOOOOE+OO 


-•3B1573E+C0 
- ♦ 199623E+00 
55*056918 
54; 982847 
2*854955 


(1 MAG) 

-♦696000E+01 
- • 194000E+ 02 
-.349O0OE+O1 

* 208000E+00 

a MAO) 

• OOOOOOE+OO 
« OOOOOOE+OO 

• OOOOOOE+OO 
♦OOOOOOE+OO 


•1707 04E+00 


UlSH ret CHANGE? 

? 1 

ENTER A NEW DESIRED VECTOR i 
? -.306 .344 -4901 -,097 4.7V .0709 1*49 
COMPLEX vd: 


-■ ♦ 306 OOOE+OO 
-■ * 970000E-01 
ACTUAL VECTOR ) 


- ♦ 98) OOOE+OO 
•70V000E-01 


♦ 479000E+01 
- • 500000E+00 


• 149000Ef 01 


** ♦ 414202E+ 00 
• 203737E+00 
ERROR VECTOR *♦ 


- « 700A04E +00 
♦ 1 2741 9E+00 


• 483360E+01 
- • 493077E +00 


♦139953E+01 


. 10G202E+00 - • 280396E+00 

- ♦ 3 0 0 7 3 7 C LOO - ♦ 5651 94E-01 

LENGTH OF THE DESIRED VECTOR ~ 
LENGTH OF THE PROJECTED VECTOR 8 
LENGTH OF THE ERROR VECTOR 
IS THE ERROR ACCEPTABLE? 


• 4360 1 9F -01 

• 692320E-02 
5*15700 7 
5*130951 

♦440622 


► 904665E-01 


EIGENVECTOR V 4) (REAL) 

- • 3 06 OOOE+OO 
901 OOOE+OO 

♦ 479000E +01 

♦ X 49000E + 01 
= B «»»K(;UNTEKr3 OF ' CURRNT ' 
MATRIX V J 


(1 MAG) 

- *344000+: +00 

• 9 70 00 OE *“01 
-* 709000E-01 

♦ 50 00 00 E +00 
DATA FILE INCLUDE* 


♦ 40 2 2 7 45410 7 7 2 9 E 1 0 2 
.350 9 X 5 3 1 5 4 4 1 0 2E f 0 0 


. 7 063 9 9 3 3 0 5 0 2 4 5 E f 0 1 


\ 1420: 


♦ 1 2 X 0 93002305 i 2E +0 X 


*Or«i 'wr-rtr 


• • 70060- 


« 6 960 OOF +01 
• 786399E+01 
- • 903993E + 00 


• 344000E+00 
.35095 SET 00 
- « 69 1532E-Q2 


3 

1 9535960 1 E tOO 
4 1 6 0 7 3 1 49 E 1 0 0 
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20373AAA88AA23Efoo 


» 5 278X5732 522C0E+02 
. 1274 5 ?;jya«s:!5SEf oo 


. 373A A34459A744E f 0 1 


.4833601 8720 A34E +01 


4 .8OOi295702A792EfOi -.83720A9A29A457E-02 

-♦493C7A797K2403E+CO 

UXSH TO IJIOPLAr THE NORMAL IZEtt EIGENVECTORS? 
r 0 

(jA.CN MATRIX + • 

.1 2 
4 

X .131 9008 A08S2A2E +0 1 - . 1 AS 1 4321 1 3A434E + 0 i 

- ♦ 3 72387224888A5E+01 


. 1 399S3354S9721 E + 0 5 


2 - • 38A30A38534AA9E+01 

.35 589348A48OA7E+02 
MATRIX AHATJ 

1 

4 


.1A927799A49577E+00 


.A0A347O4908539E+0O -.S8O853877AG71SE+01 


21 77751 3887797E+01 
230108AB232135E+01 


. 39SS22S0 1 39 A 1 AE +0 1 


74293020740597E+00 - ♦ 1 053254420881 VE + OX 

i 20 i 5542931 493E f 0 1 


3 .74AlA871b8A85AEf00 
- ♦ 849SA35048581 7E+0X 

4 . 23428 120A38957E+00 
-•43083S1955289SE+01 

WISH TO EXIT FROM THIS MODE? 


•2832904350A728E+00 


♦ 3528573055A458E+00 


18080S4074970AE+00 


-.31 A48892353 5 75E f 00 


. 5393577A4951 X AE+OO 


•841S972778599AE+00 


XX*********XXX*X*t*t*%*%%** EXITING MODE 3 #*###*****#**)|c****Cj|c***#* 

GAIN MATRIX FI 

X 2 3 4 

1 .131901E+05 - ♦ 5 651431: +01 .1A9278E+00 -.172387E+0X 


-,3U630AET03 , A0A347E+00 -.580854E401 .31X893E+02 


MATRIX VFS 

l 2 ;? 4 


i 

.1943J54E+02 

.333432E+00 

-♦ 102V0AC+00 

- • 3A0320E—01 


.A7A4A0E+01 

. 72472AE+01 

.752159E-05 

• 97A052E-01 

3 

•10A152E+01 

— 4 V74030E-01 

4 203 980E+02 

4 X X3337E+00 

4 

WANT 

- . 424*22 SE+OO . <43S* 5SK+00 

TO CONTINUE SEARCH? 

-4 8S43 SOF+Oi 

— ♦ 39 X 375E+ 0 X 

WEIGHTING' CONSTANTS 
WEIGHTS. 

X 

RO EIGENVECTOR 

2 

IMPR0VEMENT^*^***J(£***3|t*Jjc#***)tcjit* 
3 A 

X 

.i OOOOOE+OX 

• 1 OOOOOE+OX 

• X OOOOOE+OX 

• 100000E+03 

2 

•100000E+01 

• X OOOOOE+OX 

. X OOOOOE+OX 

♦3 OOOOOE+OX 

3 

.lOOOOOE+OX 

« X OOOOOE+OX 

4 X00000E+01 

4 100000E+01 

4 

WISH 

. 3 OOOOOE+OX 
TO CHANGE? 

* X OOOOOE+OX 

. X OOOOOE+OX 

4 1 OOOOOE+OX 


? :i 
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ENTER NEW VALUES } 

? x X XFIS XFZS X i XES XES If. 7 IE 7 1 X XE7 1E7 X X 

cost- ,343yy3E+oo 

GRADIENT SEARCH ROUTINE , SET SEARCH PARAMETERS} 

DEFAULT VALUES ARE} 

$ OF STEF'S > N= X STEP S I Z E r J.< ~ . XOOOOOE-OX DMIN- . I OOO OOF- 06 

WISH TO CHANGE? 

? 0 

hew cost- .;5()S:»s:>e+ob 
COST FUNCTION- . 30S252E+08 
WISH 10 CONTINUE THE SEARCH? 

? X 

GRADIENT SEARCH ROUTINE. SET SEARCH PARAMETERS} 

DEFAULT VALUES ARE 5 

* OF STEPS. N= X STEP SIZE.D- , XOOOOOE-OX DM1N-- . X00000E-06 

WISH TO CHANGE? 

? 1 


ENTER NEW VALUES} 

? XO .X .XF.-S 
NEW COST- . X284V3E+0V 

LAST 87 EP NOT ACCEPTED.! 

STEP SIZE REDUCED TO} .500000E-0X 
NEW COST- .4H'ty;WE+08 

LAST STEP NOT ACCEPTED! 

STEP SIZE REDUCED TO} .2S0000E-0X 
NEW COST* .31730SE+08 

LAST STEP NOT ACCEPTED! 

STEP SIZE REDUCED TO l . X23000E-0X 
NEW COST- . L’yXHSXE+OS 

NEW COST- .31 7305E+08 

2 STEPS WITH PRESENT GRADIENT AND DM IN- 
LAST STEP NOT ACCEPTED! 

NEW COST- .2?2885E+08 
LAST STEP NOT ACCEPTED! 

STEP SIZF1 REDUCED TO} .62S000E-02 
NEW COST = . 287M7XE+08 

NEW COST- .2V288}5F:+08 

2 STEPS WITH PRESENT GRADIENT AND DM IN- 
LAST STEP NOT ACCEPTED! 

NEW COST- .282891E408 
NEW COST- • 28S721E408 

2 STEPS WITH PRESENT GRADIENT AND DM1 N- 
l-AST STEP NOT ACCEPTED! 

NEW COST- . 277427E+08 


NEW COST- .279S87E408 

2 STEPS WITH PRESENT GRADIENT AND DMIN- 
LAST STEP NOT ACCEPTED! 

NEW COST- .272717E408 
NEW COST- . 27S343E408 

2 STEPS WITH PRESENT GRADIENT AND DM IN- 
LAST STEP NOT ACCEPTED! 

NEU COST- . 267448E+09 
NEW COST- . 27033XE408 

2 STEPS WITH PRESENT GRADIENT AND DMIN- 
LAST STEP NOT ACCEPTED! 

NEW COST- .24324XE408 
NEW COST- .244434E408 

2 STEPS WITH PRESENT GRADIENT AND DMIN- 
I- AST STEP NOT ACCEPTED! 

NEW COST- . 2S86X2E408 
NEW COST- .2M890E408 

2 STEPS WITH PRESENT GRADIENT AND DMIN- 
LAST STEP NOT ACCEPTED! 

NEU COST- .2S447VE408 
NEU COST- . 2S8223E408 

2 STEPS WITH PRESENT GRADIENT AND DMIN- 
l-AST STEP NOT ACCEPTED! 


4 X25000E-0XWERE TAKEN 


4 4 250 00 E -02WERE 7 AKEN 


. 425000E-02WERE 7 AKEN 


. 425000E-02WERE 7 AKEN 


. 425000E-02UERE TAKEN 


. 425000E-02WERE TAKEN 


. 4 2 SO 0 OE 02WERE TAKEN 
. 425000E-02WERE TAKEN 


. 425000E -02VJERE TAKEN 
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N'fc'U LUST= .2S0227Ff03 
COS T FUNO f I {Ju- • 25022 /!-:»•<:» 

UJ-SH 10 C0KTJN8JE THE SEARCH? 

"THE SEARCH IS CONTINUED AND A FINAL COST FUNCTION 

CRABXEHT SEARCH R«UTI.\*E» SET SEARCH PARAMETERS! 


IS CALCULATED."- 


C&PAU1L f V ril.MJr.f) ARK* 

* am" STEPS# N- XO STEP SIZE»to= * 625000£-02 


* .tilWHrtrtr-Oi; 


COST FUNCTION® • 1 SS 136 E +07 
MATRIX V i 


X * 48266531 19023 SE +02 

* 35093058 134550F+00 


2 , 12225367214937 E +01 
• 20373661449714 E +00 

3 . 12641703663988 E +02 
* 127431 £ 57581 04 E +00 


4 , 79780 X X SO i 0709 E+ 0 1 

- * 49307 X 80 X 4 X 918 E +00 
MATRIX VPS 

x 2 


1 *X 94545 F +02 

2 .676321 E +01 

3 * 451053 E +00 


• 344340 E +00 

* 725272 E +01 

• 33263 f;E +00 


4 - * 672359 E- 0 X « 675862 C 400 

WANT TO CONTINUE SEARCH? 

? 0 

(JAIN MATRIX 

X 

4 


• 78674405727293 K+ 0 X 
* X 980200 30 20 286 E + 0 2 
• 3 G 436583280420 E +01 
• 54 S 29 S 36274 S 72 E -01 


- « 4 X 4 X 9060399976 E +00 
70060734651333 E +00 
• 48336001377353 E+ 0 X 

• 139953453046541: +01 

« 


3 


4 


1029031*400 - * 360261 E -01 

* 7521 52 E -01 ♦ 9 76084 K -OX 


« 20 X 900 E 402 * 1 13366 C +00 

- . 8 S 4 1 7 9 E + 0 1 - * 391 3 7 6 K 4 0 X 


3 


X . X 3188024894946 K +01 

17236 1 2857 1337 E+ 0 X 

2 - * 38537988320367 E +01 

• 311 860 X 973 SO X 9 El 02 
MATRIX AHATI 

X 2 


- . X A 804 4 8 X 4 9 X X 3 7 FI I OX 


* X 6 9 3 2 8 4 4 0 8 7 4 1 1 C 1 0 0 


♦ 58332789768427 E 4 00 


- * Hfi X 0 X X 72833928 E +0 3 


3 


4 


i -. 217757 E +01 

• 395400 E+ 0 X 

-• 180862 E +00 

• 230072 E 40 J 

2 - • 742386 E +00 

-• 105432 E +01 

-• 316572 E +00 

* X 2 G 147 K 40 X 

3 * 743559 E +00 

• 289643 E +00 

•S 39798 E +00 

- * 8497751 * 40 X 

4 . 232922 E +00 

PAUSE 12 

• 356109 E +00 

• 84182 SE +00 

- * 430791 E+ OX 

F XFORMEDi 




X 

2 

3 

4 

6 

7 




j 

* 36624 if+OX 
•294204K-02 

-.3iu?0:>0E+0i 
. OOOOOOE+OO 

~ ♦ 17S6S0E+00 

-•371862E+00 

-.S371H2E+00 

*> 

•1H391XE+01 
- * X X9439K400 

-.<S034a6E+0X 

.OOOOOOE+OO 

* 2697864.40 X 

•434680E401 

- ♦ 773670E-02 



WISH TO itTSPLAY TILBA SYSTEHr 

? 0 

NS= 7 NC|j= 4 
-ALAM12T : 

i 2 3 4 

* 


♦744000E+00 

•320000E-0X 


-•337000E+00 
• X X2000E+01 


-.200000I>0X 

.OOOOOOE+OO 


•OOOOOOE+OO 
• OOOOOOE+OO 


3 .. 1S4000E+00 -.249000E+00 *9960001:400 ♦ OOOOOOEKK) 

YOU MUST SELECT MUTCH 3 OUTPUTS WILL BL USED TO FEED THE OBSERVER 
THE OUTPUTS CORRESPOND TO COLUMNS IN -AI.AMi.2T 
SELECTED COLUMNS WILL FORM A MATRIX THAT MUST Bl! NON&XNGULAK 


ENTER OUTPUTS TO BE USED (INTEGER) 
? i 2,7 


M*********^******* ASSIGN EIGENVALUES EON OBSERVER********** 
******************** MO BE 2$ EIGENVALUE ASSIGNMENT ******************** 

********** ENTER OR CHANGE EIGENVALUES* 


PREVIOUS VALUES? 

? 0 

LAMBDA ii 

REAL 23 -♦ 150000E+0X I MAG* •XSOOQOC+OX 

WISH TO CHANGE? 

? 1 

ENTER NEW VALUE ( S ) i 
? -5 0 

NEXT eigenvalue: 

PREVIOUS VALUES? 

■? a 

LAMBDA 

PEAL = - « XSOOOOE+OX I MAG 32 -♦ X15OOO0E+OX 

WISH 10 CHANGE? 

? .t 

ENTER NEW VALUE (S) i 
? -6 0 

NEXT eigenvalue: 

PREVIOUS VALUES? 

? 0 

LAMBDA .U 

REAL 52 - • 200000E+0X IMAli« « XOOOOOE+Ol 

WISH TO CHANGE? 

? 1 

ENTER NEW VALUE! S) t 
? -7 0 

WISH 10 EXIT PROM THIS MODE? 

? X 

*************************** EXITING MODE 2 ************** 


*************** ASSIGN EIGENVECTORS FOR OBSERVER************** 
******************** MODE ^EIGENVECTOR ASSIGNMENT ***************: 
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********** ENTER OR CHANGE .EIGENVECTORS* 


PREVIOUS VALUES? 

? 0 

EIGENVECTOR V XS (REAL) 

•4B7000E+02 

-.X37000E+0X 

tX2400<iE+02 

WISH TO CHANGE T 

? X 

ENTER A NEW DESIRED VECTOR S 
? 1 0 0 0 0 0 

DESIRED VECTORS 


♦XOOOOOE+OX • OOOOOOE+OO ,0000001100 
ACTUAL VECTOR S 


. XOOOOOE+OX • 00000 OE+OO .OOOOOOE+OO 

ERROR VECTOR J 


♦ OOOOOOE+OO . OOOOOOE+OO 
LENGTH OF THE DESIRED VECTOR - 
LENGTH OF THE PROJECTED VECTOR® 

LENGTH OF THE ERROR VECTOR s 
IS THE ERROR ACCEPTABLE? 

? 1 

NEXT EIGENVECTORS 
EIGENVECTOR V 21 (REAL) 

« 4G7000E+02 
- . i 'S/OOOE Un 

.124000E+02 -.349000E+OX 

WISH TO CHANGE? 

? X 

ENTER A NEW DESIRED VECTOR S 
? 0 0 X 0 0 0 

DESIRED VECTORS 

.OOOOOOE+OO ♦ 100000E+01 .OOOOOOF 100 

ACTUAL VECTOR S 

♦OOOOOOE+OO . XOOOOOE+OX .OOOOOOE+OO 

ERROR VECTOR S 

. OOOOOOE+OO , OOOOOOE+OO 

LENGTH OF THE DESIRED VECTOR « 

LENGTH OF THE PROJECTED VECTOR® 

LENGTH OF THE ERROR VECTOR 
IS THE ERROR ACCEPTABLE? 

? X 

NEXT EIGENVECTORS 
EIGENVECTOR V 3S (REAL) 

- * JOAOOOE+OO 
- ♦ 9HX000E+00 
♦4790 OOE+OX 

WISH TO CHANGE? 

? 1 

ENTER A NEW DESIRED VECTOR S 
7 0 0 0 0 % 0 

DESIRED VECTORS 


♦OOOOOOE+OO .OOOOOOE+OO .XOOOOOE+OX 
ACTUAL VECTOR S 


.OOOOOOE+OO .OOOOOOE+OO .XOOOOOE+OX 

ERROR VECTOR S 

. OOOOOOE+OO . OOOOOOE+OO . OOOOOOF+OO 

LENGTH OF THE DESIRED VECTOR ® X. 000000 

LENGTH OF THE PROJECTED VECTOR- 1.000000 

LENGTH OF THE ERROR VECTOR =■ .000000 

IS THE ERROR ACCEPTABLE? 

? X 

M «. a ...»(:DNrENTS OF 'CIJRRNT* DATA FILE INCLUDES 
MATRIX V S 

1 2 3 


.OOOOOOE+OO 
X. 000000 
1.000000 
.000000 


< IMAG ) 

♦ 344000E+00 
-.970000E-0X 
• 709000E-01 



(IMAG) 

♦ <f.9*000F.+(U 

• X 94000E+02 
• 349000E+OX 


X 0 0 0 0 0 0 0 0 0 0 00 0 E+O i 


OOOOOOOOOOOOOOE+OO 


0 00000 OOOOOOOOE+ 00 
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2 .oooooooooooooo+:+oo . ioooooooooooooe+ox 

3 . OOOOOOOOOOOOOO+3+OO ♦ 00000000000000^ f 00 

WISH TO DISPLAY THE NORMALIZED E F OEM VECTORS? 

? 0 

CAIN MATRIX Ki 

1 2 

A . i 981862327/969E+02 . 7663561 1342999E1 Ox 

2 - ♦ 36624643633 339E+00 ♦ 16743269601 H77f + 02 

;< -.320380720433 62+.+03 .300108003 6 .5 3 0 3 13 I 0 X 

MATRIX AH ATi 

1 2 

i - . 3000 000000000 XK+ 03 -.319/44 2 3 X 0 9 2 0 3 E - X 3 


.00000000000^00+3+00 
. 10000000000000+3+01 

3 

-♦ 18371670273866E+00 
. 3306 1 9 X 3060189+3 -02 
- . 6 9 V 8 0 70 6 6 8 768 2 E + 0 i 
3 

. 1 7 7 6 3 3 6 8 3 9 400 3 C - 1 4 


2 * 33527 136788005E- 14 - . 600000000000021:+ 01 - .33311 151 23125fifc> 16 

3 . 14210834 7 1 3 2 02+3-1 3 , 20421709430404E-13 ~ . 70000000000000+3+01 

WISH TO EXIT + ROM THIS MOPET 

? i * 

EXITING MODE 3 ***#*: t%%*.*X%#%***%*****%% 

? 

Li 


1 .1981868+02 - . 366246E+00 - ♦ 320389E+0.1 .OOOOOOE+OO 

2 . 766536+3+0 i . 167433E+02 .30010813+01 .OOOOOOE+OO 

3 -.185/17+3+00 *330619+3-02 699007+3+ 01 .0000008+00 

MATRIX Ei 

1 2 3 

1 -.300000+3+01 .333271+3-14 . 142109E-13 

2 -.319744+3-13 -.600000+3+01 .284217+3-13 

3 . 177636+3-14 -.333112+3-16 -.700000+3101 

MATRIX RJ 

1 2 3 4 

1 .286629+3+00 -.292894+3-01 -.283113E+02 -.371862+3+00 

2 • 1 84378+ +00 .710373+3+00 - . 171 387+3+02 .434680+3+01 

+-’AUSL 20 

■? 

MATRIX (iJ 

12 3 4 

1 -.891838+3+02 .234728E+01 -.17810813102 . J23747E+00 

2 -.421606+3+02 -.836933+3+02 .576131E+02 -. 11384 2E+00 

3 -.329284+3+01 -.34681613-01 .484813E+02 .270126E+00 

MATRIX TBJ 

1 2 

1 . 200000+3+02 . OOOOOOE+OO 

2 *000000+3+00 .250000+3+02 

3 .OOOOOOE+OO .000000+3+00 

ENTER TKitrAL CONirx T IONS'” FOR" OR I (5 1 FOTT‘T>TfTTEC 
? 0001000 

USE TH+3 E8LL8WIN8 INITIAL. CONDITIONS IN TIME RESP# 

XCO) i 


1 .00000000000000+3+00 

2 ♦ 00000000000000+3+00 

3 .00000000000000+ +00 


4 


10000000000000+3+01 
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ti ♦ ooooooooooooooe+oo 

A . ooooooooooooooe+oo 

? . ooooooooooooooe+oo 

H .OOOOOOOOOOOOOOE+OO 

v . oooo ooooooooo oe+ oo 

-10 .OOOOOOOOOOOOOOE+OO 

WISH T Cj DO TIME RESPONSE? 

? i 

MODE 4 .TIME SIMULATION %%%%%%%%%%%%%%%%%%*%*%* 

********** choose simulation options: 

-ENTER.* 1 TO SIMULATE t.’Al t 2 TO SIMULATE CAHAT.TtC! FOE l.'A'IU l)i 
? 2 

EUCER 0 TO SIMULATE OUTPUTS » 1 TO SIMULATE STATE VAR] Aft FSI 
? A 

ENTER SIMULATION TIME t (REAL. NUMBER IN SECONDS) i 
? 1 0 

ENTER NUMBER OK POINTS TO BE CAL OIJI. A T ED % ( 200 MAX ) i 
? 200 


SPECIFY THE INITIAL CONDITIONS: 

X i((J){ 

? 0 

X 2(0) : 

? 0 

x ;uo)i 
? o 

x 4 < o ) : 

? i 

X S<0) i 

? o 

x 6<o>: 

? o 

X 7(0) J 
? 0 

x o<o>: 

r 0 

x y( 0 ): 

? o 

X10(0) : 

? o 

CHOOSE INPUT OPTIONS : i FOR NO INPUT » 2 FOR A STEP INPUT r 

3 FOR A RAMP r AND 4 FOR A TRUNCATED RAMP: 

INPUT OPTION FOR U 1 : 

? 1 

INPUT OPTION FOR U 21 
? X 

ENTER 0 FOR SO DISPLAY COLUMNS » I FOR 129 COLUMNS ! 

? 0 

ENTER 0 FOR INDIVIDUAL AND 1 FOR MULTIPLE PLOTS: 

? O 

CO YOU WISH TO SET THE MIN-MAX RANGES FOR THE AXES? 

? X 

ENTER MIN X,MAX X * MIN Y t AND MAX Y VALUES ♦ 

T 0 .12 -.005 .0415 

POSITION PAPER AT TOP OF FORM AND TYPE ANY INTEGER 
YOU MAY ADD A SHORT NOTE (20 CHARACTERS. ) 

THE TIME RESPONSES ARE NOU PLOTTED. RESULTS ARE SHOUN IN FIGURES 5. lb-5. 25 
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